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1.0 SUMMARY 


1.1 Contract Description 

This contract studies unique optical computing concepts which use nonlinear 
optical phenomena to perform matrix multiplication and to provide reconfigurable opti¬ 
cal interconnection. The study focuses on the use of real-time holography in nonlinear 
media such as photorefractive crystals for optica! computing. 

1.2 Scientific Problem 

By incorporating the parallel nature of optics in nonlinear media, it is possible 
to perform parallel matrix multiplication using four-wave mixing. In addition, the 
dynamic holography in nonlinear optical media provides a natural candidate for the 
reconf igurable interconnection. The general problem in this program is to generate and 
investigate nev. concepts which use these nonlinear optical phenomena for optical 
computing. 

Specifically, this program investigates experimentally and theoretically the 
multiplication of matrices using optical four-wave mixing in nonlinear media, and the 
possibility of using such matrix multiplication and wave mixing for reconf igurable 
interconnection. 

1.3 Progress Summary 

There are several areas of significant progress achieved under this contract 
that are directly related to the development of optical matrix multiplication and recon- 
figurable optical interconnection. These include: 

1. First experimental demonstration of parallel matrix-vector multiplier 
using optical four-wave mixing in a barium titanate fBaTiO-j) crystal. 

2. First experimental demonstration of summation process inside the non¬ 
linear media in matrix-vector multiplications. 

3. Experimental demonstration of 2 x 2 matrix-matrix multiplication using 
optical four-wave mixing in a BaTiO-j crystal. 
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4. Experimental demonstration of matrix-matrix multiplication using color 
multiplexing. 

5. Experimental demonstration of matrix-matrix multiplication using convo¬ 
lution. 

6. Experimental demonstration of matrix-vector multiplication using a 

spatial light modulator and a phase conjugator. 

7. Development and experimental demonstration of a new concept of 

reconfigurable optical interconnection using photorefractive holograms. 

8. Experimental demonstration of high-efficiency interconnection using 

BaTiO^ crystals. 

9. Development of matched amplification at the Fourier plane to achieve 
maximum efficiency. 

10. Experimental demonstration of high efficiency in photorefractive optical 
interconnection using matched amplification. 

1 1. Experimental demonstration of reconfigurability using a liquid crystal TV 
in conjunction with a photorefractive barium titanate crystal. 

1.4 Publications and Presentations 

1. "Reconfigurable Interconnection Using Photorefractive Holograms," Proc. 
SPIE, Vol. 1151, Paper No. 03, to be presented at SPIE's Annual International 
Sumposium on Optical and Optoelectronic Applied Science and Engineering 
(August, 1989, San Diego, CA). 

2. "Energy Efficiency of Optical Interconnections Using Photorefractive 
Holograms," submitted to Appl. Opt. (1989). 

3. "On Image Amplification by Two-Wave Mixing in Photorefractive Crystals," 
submitted to Appl. Opt. (1989). 

4. "Photorefractive Nonlinear Optics and Optical Computing," Opt. Eng. 28(4), 
328 (1989). 

5. "Energy Efficiency of Optical Interconnection Using Photorefractive Dynamic 
Holograms," Optical Computing, 1989 Technical Digest Series, Vol. 9 (Optical 
Society of America, Washington, D.C.), pp. 128. 
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fa. "Lnergy efficient Optical Interconnection Using Dynamic Holograms in 

Photorefractive Media," OSA Annual Meeting. 1988 Technical Digest Series, 
Vol. I I (Optica! Society ot America, Washington, D.C., 1988), pp. 178. 

7. "Optical Interconnection Using Photorefractive Dynamic Holograms," Appl. 
Opt. 27, 2093 (1988). 

8. "Optical Matrix-Vector Multiplication Using a Spatial Light Modulator and a 
Phase Conjugator," Spatial Light Modulators and Applications, 1988 Technical 
Digest Series, Vol. 8 (Optical Society of America, Washington, D.C., 1988), pp. 
208. 

9. "Optica! Matrix-Matrix Multiplication Using Multicolor Four-Wave Mixing," 
Proc. SP1F, Vol. 88 I, paper 38, presented at OF/LASF'SS (Jan. 10-17, 1988, Los 
Angeles, C A). 

1C. "Optica! Matrix-Vector Multiplication via Four-Wave Mixing in Photorefrac¬ 

tive Media," Opt. Lett. J 2, 138 (198 7); Opt. Lett. 12 , 37 3 (1987). 

11. "Optical Matrix-Vector Multiplication Using Four-Wave Mixing," Paper MM5, 

OSA Annual Meeting (Oct. 1986, Seattle, W A), J. Opt. hoc. Am. A3 (13), lb 
(1986). 

1.5 Pa te n t Disclocires 

"Multicolor Matrix-Matrix Multiplier with Parallelism," M. Khoshnevisut., A.L.T. 

Chiou and F. Yeh, Rockwell Patent Disclosure 87SC37. 

"Matrix-Vector Multiplier using Photorefractive Phase Conjugators," P. Yeh, A.L.T. 

Chiou and M. Khoshnevisan, Rockwell Patent Disclosure 87SC96. 

"Optical Matrix-Matrix Multiplier with Parallelism by Spatial Convolution via Four- 

Wave Mixing," A.L.T. Chiou, Rockwell Patent Disclosure 87SCfa3. 

"Reconfigurabie Optica! Interconnect using Dynamic Holograms," P. Yeh, Rockwell 

Patent Disclosure 87SC54. 
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2.0 TECHNICAL DISCUSSION 


The technical problem addressed in this contract is the study of optical matrix 

multipncaMon and recordigurable optical interconnection. By incorporating the parallel 

I 2 

nature o! optics m nonlinear media such as photorefractive crystals, ’ it is possible to 

perform parallel matrix multiplication using four-wave mixing. In addition, the dynamic 

holography in nonlinear optical media provides a natural candidate for the reconfiguratle 

3 L, 

optical interconnec 11 on. 

Specifically, mis program is aimed at the utilization of nonlinear optica! 
phenomena such as phase conjugation,'* ' two-wave mixing,and tour-wave mixing^"' 7 
m nor,linear media ’ to perform trie m; trix-vector and matrix-matrix multiplications 
afid reconf gurarable optica! interconnection. Photorefrartive crystals such as BaTiOj 
and strontlum o.-.rium niohate (SBN) we-e used as the nonlinear media. 

c«ei.era.'!y speaking, tour-wave mixing is a nonlinear optical process in wfuch 

three input waves mix to y leid a fourth, wave. In phase-matched four-wave mixing, the 

three input waves consist of two < ounterpropagating pump waves, E j and Tp- and an arbi- 

(hi 

:rar> prose wave, L -M: three couple through the third-order susceptibility, > , to 

yield a ! 0 ' ; ’ w.-.ve. L 4> w‘ n is proportional to the produi t of IS j, 11 2 «*nd complex coiiju- 
. . , S 

cl ■ Oi I 7 , d f i*_ < pi"' JL* V» ' . \ f t»l. cl S 


Such a tour wave 1mx11.it c at be understood in terms of the recording and readout proces¬ 
ses whii h. occur in hologr ap f, y. In UK- nonlinear media, one of the pump waves and the 
probe wave form an interlerenc e pattern which, in turn, induces an index grating. The 
other pump wave is diffracted from this grating and generates the fourth wave. The 
formation of the grading (or hologram) and readout processes take place at the same 
time. Thus, four wave mixing is sometimes referred to as a real-time holography. 

Bv using four wave mixing in nonlinear media, multiplication ot signals can 
take place m a subpic osecoi.d tune scale. In addition, it we use the pat a I lei nature of 
optit a. wa..',, each wave 1 an i arrv suatial informatioii for tfic* purpose o' image* process¬ 
ing u-.ing the two tram.ver -e dunm.ioru. Sue! an applu at ion of four - wave mixing to 
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real-time imag processing was recently reported using BSO crystals.^ In what follows, 
we will describe some unique concepts which use the two transverse dimenoions to carry 
the matrix information for the purpose of matrix multiplication and reconfigurable 
optical interconnection. 

2.1 Ma trix Matrix Multiplication 

The matrix multiplication^ ^ between two N x N matrices can be stated as 

fol lows: 


C = AB 


(I) 


where 


\ • , U, • 

IX k J 


( 2 ) 


Note that a matrix multiplication consists of two mam operations, a parallel multiplica¬ 
tion and an addition. 

Referring to Fig. ia, let us consider a four-wave mixing configuration which is 

1 3 

suitable for matrix multiplication. Beam 1 and Beam 2 contain the information about 
the two matrices A(x,z) and B(>:,y), respectively. Beam 1 is propagating along the y-axis, 
and Beam 2 is propagating along the x-axis. These directions are chosen for the sake of 
clarity in introducing the concept. They are not the only directions for matrix 
multiplication. Also, these matrices can be either continuous or discrete. In the discrete 
case, each beam consists of a matrix of beamlets, as shown in Fig. lb. 

In the nonlinear medium, these two matrix-carrying beams form an interfer¬ 
ence pattern. As a result of the nonlinear response of the medium, a volume grating is 
formed. This grating contains information about the product of the matrix elements of 
these two matrices, and can be written as: 


r.n = 


n ? A(x,z)B*(z,y)e 




+ c.c, 


(3) 
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Clx yl - Al* »)B* li.yld/ 


*C«t 3 1 640 

READOUT 

BEAM 


(b) 



Fig. I Matrix-matrix multiplication using four-wave mixing. 

- j?.s the difference of the wavevectors of the matrix-carrying beams and cc repre- 
complex conjugation. The parameter r >2 is the Kerr coefficient and is proportional 
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to the third-order susceptibility of the medium. Note that the nonlinear response of 
the medium serves the purpose of parallel multiplication. 

The volume grating is then read out by a third beam which can simply be a 
plane wave. The diffracted beam consists of the integrated contribution from each part 
of the grating along the beam path, and thus can be written 

C(x,y) <* J A(x,z)B*(z,y)dz (4) 

where the integration is carried out along the beam path. Note that the integration com¬ 
pletes the operation of matrix multiplication. The information about the product of 
these two matrices is novn, impressed on the transverse spatial distribution of the dif¬ 
fracted beam. 

Due to the phase matching requirement, the readout beam must be incident 
along the directions which satisfy the Bragg diffraction condition to achieve high effi¬ 
ciency. In anisotropic nonlinear media, the polarization states, as well as the direction 
of propagation, can be chosen such that the largest of the nonlinear susceptibilities is 
fully utilized. 

The four-wave mixing can be either degenerate or nondegenerate. In the 
degenerate case, all the beams have the same frequency. In the nondegenerate case, the 
frequencies of the beams can be slightly different. This may be useful for the purpose of 
separating the diffracted beam from the undiffracted portion. 

To illustrate the information capacity of such a nonlinear optical matrix 
multiplication, let us consider a four-wave mixing process using an Ar ion laser at 48S0A 
in a medium of a I cm cube. The grating space is of the order of 0.5 uni. Thus, 10 um x 
10 urn is enough for each pixel of information. In other words, a 1 cm cube is cable of 
handling 1000 x 1000 matrices. \X ith a material response time of 1 ns, such a matrix 
multiplication has a potential data throughput rate of quadrillion bits per second 
(10 15 bits/s): 
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2.2 Matrix-Vector Multiplication 

Another scheme as shown in Fig. 2 is suitable for matrix-vector multiplica- 

11 - 

tion. Here, as an example, let us consider a discrete case in which we need to carry 

out the multiplication of an N-element vector and an N x N matrix. The vector is fanned 
out into N-rows of identical vectors. These N x N small beams are directed to a non¬ 
linear medium. The matrix which also contains N x N small beams is also directed to the 
medium in such a way that each beam of the matrix is counterpropagating in a direction 
relative to the corresponding beam of the vector. Thus, in the medium, there are N x N 
spatially separated regions which are pumped each by a pair of counterpropagating 
beams. Now, N x N probing beams are directed into the medium in such a way that each 
probe beam will propagate through an intersection region. The probe beams will be 
"plane wave" beamlets propagating in parallel. As a result of the four-wave mixing, each 
probe beam will generate a phase-conjugated beam which, within a proportional factor, 
can be written M(i,j)a(j). By using a cylindrical lens, a summation over j can be ob¬ 
tained. Thus, we have 

b(i) = V M(i, j)a(j) (5) 

J 

where a(j) is the j-th element of the vector a and M(i,j) is the matrix element. Such a 
scheme for matrix-vector multiplication can also be used for matrix-matrix multiplica¬ 
tion by decomposing a matrix into column vectors and then multiplying the matrix with 

each of the column vectors. Using color-multiplexing, one can perform matrix-matrix 

3 15 

multiplication with N -parallelism. 

2.3 Reconfigurable Optical Interconnection 

Optical interconnection will play a key role in both the optical computing and 
VLSI systems.There are many advantages of optical interconnect. These include 
high space-bandwidth and time-bandwidth products, thus, many independent channels 
could be exploited for demanding computations. Optical processors are inherently two- 
dimensional and parallel. Optical signals can propagate through each other in separate 
channels with essentially no interaction. Optical signals can also propagate in parallel 
channels without any interference and crosstalk. In the VLSI systems, optical inter 
connect can be used to solve the problem of communication, as well as the clock 
distribut ion. 
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Fig. 2 Matrix-vector multiplication via four-wave mixing. 

Computer-generated holograms (CGH) can be used for fixed optical intercon¬ 
nect. The most general CGH consists of a two-dimensional array of subholograms. Each 
subhologram is capable of diffracting a gate output to any gate or combination of gate 
inputs. Such a fixed interconnect pattern is adequate for many applications, including 
the point operation, matrix operation, the clock distribution and a globally synchronous 
systolic processor. However, there are many situations in which a dynamic optical inter¬ 
connect is required. 

In Fourier transform and sorting, every element of an output array is affected 
by all the elements of the input array, and conversely, each element of the input array 
affects all elements of the output array. In addition, the interconnections change at dif¬ 
ferent stages of the computation. Thus, these operations require global and dynamic 
interconnection between different elements of the input array. 

In image restoration, pattern recognition, and neural network systems, the 
interconnect patterns could be data-dependent, making it impossible to foresee the 
interconnect requirements at different stages of processing without having foreknow¬ 
ledge of the input. The computational throughput of parallel processor implementing 
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these types of operation will be critically affected by the availability of a dynamic and 
global interconnect network. 

In Fourier transform and sorting, the requirement for dynamic interconnect 
can be avoided by resorting to a fixed but global interconnect pattern known as the "per¬ 
fect shuffle". In the area of image restoration and pattern recognition, there is a mini¬ 
mum amount of regularity and structure. In addition, there is a possible data and time 
dependency in the interconnect requirements. A global and reconfigurable interconnect 
network will be vital to achieving high throughput and high efficiency with parallel proc¬ 
essors implementing these operations. 

The most general interconnect system is one in which any gate output can be 
connected to the input of any gate or combination of gates (see Fig. 3). The effect of 
such an interconnect can be represented by the matrix equation 

0 = MI (6) 

w'here I is a vector representing the two-dimensional input array, M is the matrix repre¬ 
senting the interconnect, and O is a vector representing the output array. In digital opti- 
17 18 

cal computing, ’ the input array is actually the gate output array. Each matrix- 
element (i,j) is nonzero if, and only if, there is a connection between pixel j of the input 


•cs*-m»4 

NONLINEAR MEDIUM 



Fig. 3 Reconfigurable optical interconnection. 
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array and pixel i of the output array. The matrix-vector multiplication technique 
described in Section 2.2 can be used as an optical interconnect, which provides both local 
and global communication between the gate input and gate output. In the VLSI systems, 
the vectors 1 may consist of N laser beams, each containing a stream of data from a 
processor, and the vector O represents the output laser beams, each feeding into a proc¬ 
essor. The main advantage of using such nonlinear matrix-vector multiplication for 
interconnection is that the matrix M can be easily changed for the purpose of reconfigur¬ 
ing the interconnect. 

Equation (6) can be generalized such that both the input array I and the output 
array O are two-dimensional matrices. In that case, the most general interconnect 
matrix M must be four-dimensional (i.e., tensor of rank four). Each matrix element (ijkl) 
represents the connection between pixel (ij) of the output array and pixel (kl) of the input 
array. Because of the two-dimensional nature of optical waves, it is desirable to repre¬ 
sent the tensor M by a two-dimensional array of two-dimensional submatrices. Thus, the 
matrix-matrix multiplication described in Section 2.1 can be used for reconfigurable 
optical interconnect which provides the most general interconnection pattern. While this 
interconnection scheme allows complete generality, a price is paid in terms of the space- 
bandwidth product requirements on the nonlinear media. 
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3.0 PROGRESS 

3.1 Matrix-Vector Multiplication 

13 19 

In our earlier demonstration of optical matrix-vector multiplication ’ (see 
Appendix 5.1), we have used optical phase conjugation via nonlinear four-wave mixing in 
a photorefractive- BaTiO^ crystal to perform parallel multiplication. The summation 
required to obtain matrix-vector products is performed subsequently by using a c,lin- 
drical lens external to the nonlinear medium. With a slight modification of the experi¬ 
mental geometry, we have also successfully demonstrated that both the parallel 
multiplication and the summation can be done inside the nonlinear medium to achieve the 
desired matrix-vector multiplication without an external cylindrical lens. 

The top view of the new experimental geometry is illustrated schematically in 
Fig. 4. The readout beam in this case consists of a vertical column of beamlets with 
equal intensity. This beam enters the crystal at an oblique angle such that each beamlet 
traverses, inside the crystal, all the counterpropagating pumping beamlets at the same 
horizontal plane. As the proper elements of the matrix and the vector are encoded in the 
two counterpropagating beams via the appropriate masks, the phase-conjugate output of 
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Fig. 4 Optical matrix-vector multiplication using photorefractive four-wave mixing. 
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each beamlet consists of the sum of the product resulting from each individual encoun¬ 
ter. From a pure geometrical point of view, the probe beam (i.e., the readout beam) can 
conveniently be injected (from c-face of the crystal) perpendicular to the counterpropa- 
gating pumping beams. 

The experimental results for the multiplication of a 4 * 4 binary matrix and a 
4 x 1 binary vector are shown in Fig. 5. The magnified images of the masks used to 
encode the matrix, the vector and the probe are shown in Figs. 5a, b and c, respec¬ 
tively. The experimental results representing the final matrix-vector product is shown in 
Fig. 5d. The relative intensity of each output spot representing the element of the 

product vector is shown in Fig. 6. Due to the slow response time (or the order of second 

2 

for optical intensity of the order of a few tens of milliwatt/cm ) and the consequential 
sensitivity to environmental changes of the photorefractive process, the output intensity 
fluctuates significantly. - 



(a) 

INPUT MATRIX 



(b) 

INPUT VECTOR 
(EXPANDED INTO 
MATRIX FORM) 



(c) 

READ OUT 
(PROBE) 
VECTOR 



(d) 

OUTPUT 
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Fig. 5 , xperimenta! results of optical matrix-vector multiplication using phoiore- 

fractive four-wave mixing. 
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Fig. 6 Output of the detector monitoring the intensity of the optical beam represent¬ 
ing the product vector. 
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Fig. 7 Decomposition of matrix-matrix multiplication into matrix-vector multipli¬ 
cation. 

3.2 Matrix-Matrix Multiplication via Color Multiplexing 

Vte have also extended our experimental demonstration on matrix-vector 
multiplication to matrix-matrix multiplication by color multiplexing (see 

Appendix 5.2). The basic idea is to decompose the problem into matrix-vector 
multiplications, as shovxn in Fig. 7, and carry out all the parallel matrix-vector 
multiplications simultaneously by using color multiplexing. Note that each of the 
operations shown on the right-hand side of the equation in Fig. 7 is in fact a matrix- 
vector multiplication. The basic principle of color multiplexing used to encode the 
component vectors with different colors is illustrated in Fig. 8 for the case of k » 4 
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Fig. 8 Optica] matrix-matrix multiplication via color multiplexing. 

matrices. Each row vector of the matrix M| is illuminated by one color, and all the color 
components are then combined by the prism (angular multiplexing) into a single row prior 
to further expansion by anamorphic optics (not shown in the figure) to match the mask 
representing the second matrix. After proper element-by-element multiplications and 
summation (summing optics omitted) into the column, the resulting multicolor output is 
then demultiplexed into different color components that together represent the final 
product. 

Using an argon ion laser that oscillates simultaneously at five colors in the 
blue-green and a BaTiO-j crystal, we have demonstrated the principle described above for 
the case of 2 * 2 matrices. Details are given in Appendix 5.2. 

3.3 Matrix-Matrix Multiplication using Spatial Convolution 

In addition to the matrix-matrix multiplication using color-multiplexed four- 
wave mixing described in the previous section, we have also successfully demonstrated 
another scheme for matrix-matrix multiplication using spatial convolution via four-wave 
mixing. From the experimental point of view, the key difference between this approach 
and the others described above is that the nonlinear crystal is now located at the com¬ 
mon Fourier plane, rather than the common image plane, of the input matrix masks. The 
encoding scheme, as explained in Appendix 5.3, is also different from those used in the 
earlier approaches. From the conceptual point of view, matrix-matrix multiplication in 
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full parallelism is achieved by space-multiplexing via spatial convolution using 
degenerate four-wave mixing. 

3.4 Matrix-Vector Multiplication Using a Spatial Light Modulator and a Phase 

Conjugator 

We have developed and demonstrated another new scheme for optical matrix- 

vector multiplication that uses a phase conjugator (with a finite storage time) in 

conjunction with a spatial light modulator (SLM) (see Appendix 5.4). The optical system 

involved is relatively simple compared with the other approaches. Without any modified- 

2 

tion, such a scheme can also perform matrix-matrix multiplication with N /2 parallelism. 

Referring to Fig. 9, we use a SLM to impress the matrix and vector informa¬ 
tion in sequence to an input laser beam. This beam is directed toward a phase conjugator 
which has a finite storage time (a photorefractive barium titanate, for example). A 
cylindrical lens is inserted in the phase conjugate beam path to perform the summation. 


COLLIMATED 
LASER BEAM 

0 



VECTOR X 

V = A>f 

Fig. 9 A schematic diagram illustrating the basic concept of optical 
matrix-vector multiplication using a photorefractive phase 
conjugator in conjunction with a spatial light modulator. 

The principle of operation is as follows: the SLM first impresses the matrix 
information onto the input laser beam. This beam is then incident into a phase 
conjugator which stores the matrix information after a finite grating formation time. 
When the matrix information is removed from the SLM, say by turning all the pixels into 
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maximum transmission condition, the phase conjugate beam which contains the recon¬ 
struction of the matrix information exists for a finite duration. This finite storage time 
depends on the strength of the input (read) beam. During this time, if the next frame of 
the SLM carries the vector information, parallel multiplication is performed as the 
phase conjugate beam propagates back through the SLM. Here, the vector is represetned 
as a two-dimensional array of N identical column vectors, where N is the dimension of 
the vector. A cylindrical lens in the output port is used to performed the summation. 
The dark storage time during which the matrix information can be retrieved is 

determined by the photorefractive material and the pumping configuration. It ranges 

21 

from seconds to microseconds. 

The system can also perform matrix-matrix multiplication by time multi¬ 
plexing. In this case, each column vector V| (i = 1 to N) which constitutes the second 
matrix M 2 is sequentially impressed onto the beam to multipl> with the first matrix Mj 
according to the matrix vector multiplication scheme described below. To avoid the 
degradation of the information of Mj stored in the photorefractive hologram during th 
readout, it is necessary to refresh the holographic memory with M| to restore its 
diffraction efficiency. This can be done by re-impressing Mj onto the beam after each 
readout cycle. Consequently, a total of 2N clock cycles consisting of N cycles of write 
and N cycles of read, will be required to carry out the multiplication of two N ■ N 
matrices. 

Using a photorefractive barium titanate crystal as a phase conjugator in 
conjunction with a US ■ L8 magneto-optic spatial light modulator (SIGHT-MOD SMD4SI 
from Semetex Corp.), we have demonstrated the basic principle described above. The 
virtue of optical phase-conjugator ensures that precise pixel by pixel alignment is 
achieved automatically. 

3.5 Optical Interconnection using Dynamic Photorefractive Holograms 

3.5.1 P rincipl e of Operation 

During the first year of the program, we have developed a new concept of 

22 

reconf igurable optical interconnection using photorefractive holograms (see 
Appendix 5.5). This new scheme provides an optical interconnection between an array of 
laser sources and on array of detectors with a very high energy efficiency. 
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Reconfigurable interconnection linking laser arrays and detector arrays plays a 
key role in optical computing. Conceptually, such interconnection can be achieved by 
using an optical matrix-vector multiplication. 

v' = M v (7) 

where v is the input vector representing the signals carried by the array of lasers, and v' 
is the output vector representing the signals carried by the array of detectors. The 
matrix M represents the interconnection pattern. When a transparency or spatial light 
modulator (SLM) is used as the interconnection pattern, a large fraction of energy is 
absorbed by the transparency or SLM. This energy loss increases as the dimension of the 
array increases. Often, the fractional energy loss can be as large as (N-D/N, where N is 
the dimension of the array. For a 1000 > 1000 crossbar switch, the loss due to fanout can 
be as big as 99.9%. This is not acceptable in high-speed computing because signals are 
passing through the SLM billion times per second ano the energy loss can be enormous. 

Referring to fugs. 10 and 11, we describe a new method of reconfigurable opti¬ 
cal interconnection which uses the nonreciprocal energy transfer in photorefractive two- 
wave mixing to improve the energy efficiency. Figure 10 describes a one-dimensional 
case for the sake of clarity in explaining the concept. A small fraction of a laser beam is 
coupled out of the beam by using a beam splitter. This small fraction (called probe 
beam) is then expanded by using a cylindrical lens and passes through the SLM. In the 
example shown, the laser beam is to connect to Detectors b and d as prescribed by the 
SLM. The transmitted beam is then recombined with the main beam inside a photore¬ 
fractive crystal. As a result of nonreciprocal energy coupling, almost all the energy in 
the main beam is transferred to the probe beam which carries the interconnection pat¬ 
tern. The results is an optical interconnection wuth high energy efficiency. F'igure II 
describes the reconfigurable interconnection for laser arrays and detector arrays. In the 
example (a 4 . 9 interconnection) shown, Laser I is to connect to Detectors b and c, 
Laser 2 is to connect to Detectors a and d, Laser 3 is to connect to Detectors c and d, 
Laser 4 is to connect to Detector a and c. A cylindrical lens is used to focus the two- 
dimensional array of beams into a vector (one-dimensional array). As a result, Detector 
a receives signals from Lasers 2 and 4, Detector b receives signals from Laser I, Detec¬ 
tor c recenes signals from Lasers I, 3 and 4, and Detector d receives signals from Lasers 
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Fig. 10 


A schematic drawing of a 1 - N optical interconnection using dynamic 
pnotorefractive holograms. N - 4. 



Fig. 1 1 A schematic drawing of a N * N optical interconnection using dynamic 
pho tore tractive holograms. N = 9. 

2 and 3. Such a concept can be extended to interconnect lasers with detectors 
where Nj_ and Nj> are two large numbei s. 

The interconnection can be reconfigured by using a different SLM pattern. 
The energy Iosj due to SLM is no more than the reflectivity of the beam sphtter plus the 
material bulk absorption due to the photorefraclive crystals. The beam splitter can be 
chosen such that the reflectivity is small (e.g., 5% or less) so that such energy loss is 
mirnmi/ed. The hologram formation in the pnotorefractive medium will limit the recon- 
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figuration time. Once the interconnection pattern is formed inside the photorefractive 

crystal as a hologram, such a scheme is capable of providing the interconnection for high 

data rate transmission. In such an interconnection, the output of each laser can be input 

to any one or all of the detectors. Optical phase conjugation can also be used in 

conjunction with the two-wave mixing to correct for any phase aberration that may be 

23 

caused by the crystal imperfection. 

3.5.2 E xperimental Demonstration 

An experimental configuration for a 1 * 32 interconnection using this 
holographic approach is shown in Fig. 12. A Ronchi-Ruling (5C line-pairs per inch) is used 
as a 32-element binary vector mask. 5% of the collimated input beam transmitted 
through the beam splitter BS is expanded vertically by a pair of cylindrical lenses (Lcl 
arid Lc2) to illuminate the vector mask. 95% of the input energy reflected by the beam 
splitter RS serves as the pump beam which interacts with the spatially modulated signal 
beam inside the barium titanate crystal. The crystal and the mask are placed at the 
bdck-focul plane and the front-focal plane, respectively, of a Fourier transform 
cylindrical lens (Lc3). The crystal is oriented so that energy is transferred from the 
pump beam to the signal beam. The amplified signal output carrying the interconnection 
pattern is roimaged by the cy lindrical lens (Lc4) on the screen SC. 
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F:g. 12 A schematic diagram of the experimental configuration lor a 
1 - 32 interconnection using a photorefractive barium titanate 
crystal. 
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The intensity distribution at the output image plane is shown in Fig. 13. The 
experimental results is shown in Fig. 14 for the case when only three out of the 32 
channels are "ON". The table in the right hand side of the figure compares the energy 
efficiencies achieved with and without the photorefractive holograms. For the latter 
case, the beam splitter BS is removed (see Fig. 12) so that all the input energy is directed 
into the signal channel. Note that the energy efficiency is improved by more than a 
factor of 7. 


32 ELEMENT VECTOR MASK 



(bl OUTPUT THROUGH B8T1O3 CRYSTAL WITH PUMP OFF 
(cl OUTPUT THROUGH BbTi 0 3 CRYSTAL WITH PUMP ON 
(1 e . AMPLIFIED VIA TWO WAVE MIXING) 


Fig. 13 Intensity distribution at the output image plane. 


3.5.3 Demonstration of High Data Rate Transmission 

Ue have experimentally demonstrated that the holographic interconnection 

22 

using photorefractive holograms can accept very high data rate signals [see 

Appendix 5.5J even though the reconfiguration time is limited by the photorefractive 

21 

grating formation, which is typically of the order of milliseconds at modest intensity. 

To demonstrate this fact, temporal modulation was impressed on a laser beam (argon, 

X - 514.5 mu) using an acoustooptic device to simulate a signal, which is to be 
interconnected with some output. The signal used was a pulse tram of frequency 
f - 0.S33 MHz with each pulse being ~ 0.2 ps wide. This rate is clearly much higher than 
the reciprocal of the photorefractive response time. The modulated laser beam was then 
split into two beams and mixed in the crystal as described before, and the amplified 
probe beam was monitored with a piiotodetecior. The upper oscilloscope trace in Fig. 15 
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Fig. 14 Intensity distribution at the output image plane when all but 
three of the signal channels are "ON." The table on the right- 
hand side compares the energy efficiencies achieved using the 
conventional approach and the holographic approach. 
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Fig. 15 (a) Signal carried by the probe beam and (b) amplified signal carried by the 

probe beam after undergoing two-wave mixing. 


shows the input probe signal, and the lower one shows the amplified probe signal. The 
results show a steady-state response in which the temporally modulated pump and probe 
beams interact simply by diffracting off the stationary index grating that is created in 
the crystal after the photorefractive response time. As this experiment was performed 
merely to demonstrate the high signal bandwidth of the system, optimization of the 
parameters was not done, and the results in Fig. 15 can clearly be improved. 
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3.5.4 Energy Efficiency of Photorefractive Interconnection; Theoretical Calculation 
and Experimental Measurement [See Appendix 5.5, 5.6, and 5.7] 

Referring to Fig. 16, we can estimate the energy efficiency of photorefractive 
crossbar switch as follows. Let 1 be the input beam intensity and R be the reflectance 
of the beam splitter. It is legitimate to assume that the beam splitter is practically 
lossless. U e may also assume that the surface of the photorefractive crystal is 
antireflection-coated so that the Fresnel reflection loss can be neglected. Under these 
conditions, the two beams that arrive at the photorefractive crystal have energies 
I 0 < 1 - R) and 1 Rt/N, respectively, where t is the intensity transmittance of the "on” cells 
of the mask. Inside the crystal, these two beams undergo photorefractive coupling. As a 
result, most of the energy of the pump beam I (1 - R) is transferred to the probe beam 
I Q Rt/N, which contains the inteiconnection pattern. The energy efficiency can be easily 
derived and is given by" 


out 
1 i n 


tR _ 1 + m 

N 1 + m exp(-’fL) 


exp (-aL) 


( 8 ) 


where m is the beam, intensity ratio, 


m 


(i' m 

tR 


(9) 


and L is the interaction lengrh, , is the coupling constant, and a is the bulk absorption 
coefficient. For photorefractive crystals such as BaTiO^ and SBN, the coupling constant 
is very large <i.e., yL >> 1). The efficiency can be written approximately 


r; 



(1 - R) i 


exp(-cl) 


( 10 ) 


23 

C 101 78TD/b je 














Rockwell International 

Science Center 


SC5502.FR 



Fig. 16 A schematic diagram snowing the intensities of various beams in a N » N 
permutational interconnection. 

We notice that for large N, the energy efficiency is limited by the crystal bulk absorption 
exp(-aL) and is maximized by using a beam spliter with a very small reflectance R (i.e., 
R ~ 0). 

In Fig. 17, the maximum energy efficiency (with optimum beam splitting ratio) 
predicted by Eq. (8) (assuming that the insertion loss of the mask is negligible, i.e., t = 1 
and that the absorption coefficient a = 0) is plotted against the dimension of the array N 
for different coupling strength yL. The diagonal straight line represents the intrinsic 
fan-out loss. For cases where a t 0, the curve remains valid provided that the vertical 
scale is multiplied by exp(-aL). 

Using an experimental configuration illustrated schematically in Fig. 18, we 

have measured the energy efficiency of the photorefractive interconnection holograms in 

13 

a barium titanate crystal as a function of the fan-out loss. Instead of expanding the 
signal beam through a vector mask as described earlier (see Fig. 12), a neutral density 
filter (NDF) in the signal arm is used to simulate the fan-out loss. The experimental 
results are shown in Fig. 19. The solid line is a theoretical fit using Eq. (8) with aL and 
yL as the fitting parameter. The results (aL = 9.2, yL = 1.1) agree with those obtained 
from other independent measurements. 

3.5.5 Maximizing Energy Efficiency by Matched Amplification via Fourier 
Transformation (see Appendix 5. 7) 

In the experiments discussed above, both signal and pump beams are Gaussian 
beams with no spatial intensity pattern. In the interconnection applications, these beams 
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ENERGY EFFICIENCY vs ARRAY DIMENSION 



Fig. 17 The theoretical energy efficiency as a function of the dimension of the array 
as predicted by Eq. (8) for t = 1, a = 0, and an optimum R ~ 0. 


SC44329 



Fig. 18 A schematic diagram of the experimental configuration for measurement of 
the energy efficiency of photorefractive interconnection holograms in a 
barium titanate crystal. 
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N 

Fig. 19 Energy efficiency (n) of photorefractive two-bearr, coupling in a barium 
titanate sample as a function of the transmittance of a neutral density filter 
placed in the signal input arm. The transmittance is labeled 1/N to relate it to 
the fanout loss of a N > N permutation crossbar network. 

are spatially modulated. The energy efficiency of the photorefractive interconnection 
depends on the spatial overlap of the beams. We have developed the utilization of 
Fourier transform to achieve maximum beam overlap and thus to achieve maximum 
energy efficiency. 

It is important to note that although the beamlets in the signal beam and the 
pump beam are intrinsically different because of the interconnection pattern, the 
individual pixels can have identical shape (e.g., square). Complete overlap is possible in 
the Fourier domain, provided all the pixels are identical. Let s(x,y) be the aperture of an 
individual pixel, and o(u,v) be the Fourier transform. Fourier transformation of a beam- 
let of the pump beam and the corresponding column of beamlets of the signal beam can 
be performed by a lens. The resulting amplitude distributions at the rear focal plane are 
given by 


Pump: P(u,v) = 

o(u,v) 

(11) 

Signal: S(u,v) = 

!. o(u,v) exp | U n (u,v) | 

(12) 


n 
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respectively, where u,v are coordinates in the Fourier plane. The summation in Eq. (12) 
is over all the windows open in the column, and 4 > n (u,v) is a phase which depends on the 
position of the window. 


W e note that the Fourier transform of a column of the signal beam consists of 
a linear superposition of identical panerns each with a different phase factor. Such 
phase factors (exp(: yn ) are due to the relative positions of the windows in each column. 
Although these phase factors may lead to interference structures, all the energies are 
confined under the same envelop c(u,v). A.s a result of the shift invariance in Fourm* - 
transformation, each of the signal beamlets overlaps completely with the pump beam at 
the Fourier plane. Thus, maximum energy coupling is achieved. 


The experimental configuration is illustrated in Fig. 20. The output from an 
argon laser (51 ^.5 nm) is spatial-filtered and expanded to form a collimated beam with 
2 c:n diameter. A variable beam splitter consisting of a half-wave plate and a polarizing 
beam splitter cube is used to split the incoming beam into two, a pump beam and a signal 
beam. Another half-wave plate is used to rotate the polarization of the reflected (signal) 
beam by 90° from the direction perpendicular to the plane of the paper to the in-plane 
direction. A polarizer placed downstream further filters out the residual perpendicular 
component. A 10 « 10 binary matrix mask is used to encode a spatial pattern onto the 
signal beam. 

To maximize the spatial overlap of the pump and the signal beams inside the 
crystal via the properties of Fourier transform, a pump mask with a single aperture 
identical to the unit cell of the signal mask is used. The two masks are placed at the 
front focal planes of two Fourier transform lenses FLj and FL 2 (of identical focal length 
f = 25 cm), respectively. The crystal is located at the common Fourier plane of the two 
masks (Fig. 20). The shift-invariant property of Fourier transformation ensures that, 
apart from a phase factor, the diffraction pattern from each signal channel overlaps with 
that from the pump inside the crystal. The intensity patterns of the pump and the signal 
at the image plane and the Fourier plane are shown in Fig. 21. Note that the two 
intensity distributions at the Fourier plane differ significantly in fine structures due to 
multiple beam interference among the various signal channels. 
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Fig. 20 An experimental configuration for a 1-to-N * N (for N = 10) broadcasting 
network using photorefractive holograms at the Fourier domain. 
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Fig. 21 The intensity patterns of the masks for the probe and the pump beams at the 
image plane and the Fourier plane. 


28 

C101 78TD/bje 

































Rockwell International 

Science Center 


SC5502.FR 

To investigate how the energy efficiency depends on the number of signal 
channels, w'e used a variable rectangular aperture in front of the signal mask to vary the 
number of channels from 100 to 1. In each case, we also adjusted the laser power and the 
variable beam splitter so that both the pump power and the total signal input power are 
fixed at 600 u\& and 6 uW, respectively. Figure 22 is a plot of energy efficiency vs 
number of signal channels. The energy efficiency turns out to be fairly insensitive to the 
number of signal channels. The crystal orientation used in this particular experiment is 
far from optimum, and one can expect significant improvement in energy efficiency once 
optimized. 
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Fig. 22 Energy efficiency (r,) as a function of number of signal channels (N) in a 1-to- 
N > N broadcasting configuration using photorefractive holograms in a barium 
titanate sample. 

3.5.6 Contrast Enhancement by Double Passage (see Appendix 5.9) 

All of our experimental results reported earlier are obtained with photographic 
masks replacing the SLM. To demonstrate the reconfigurability and to study the 
reconfiguration time, one can no longer avoid the use of SLM which, in general, have 
relatively poor contrast ratio (of the order 10 to 100). In the parallel matrix-vector 
multiplication approach, the maximum number of fan-out channels of the optical cross¬ 
bar can be limited b> the finite contrast of the SLM. 
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As a simple example, consider a M * M permutation matrix mask with contrast 
ratio C = T j/Tq, where T j and Tq are the intensity transmittance of the "ON" state and 
the "OFF" state, respectively. We assume that the signal from each source is intensity 
modulated between zero and peak value P, the signal level received by each detector 
through the "ON" cell is PTj, whereas the noise level received through the (M-l) "OFF" 
cell is P(M-1 )Tq. The signal-to-noise ratio (S/N) is therefore given by 
S/N = T 1 /(M-1)T 0 = C/(M-1) = C/M for M > 1 . The condition S/N > 1 is satisfied 

provided that the number of fan-out channels (M) is less than the contrast ratio (C). 

A simple method to improve the contrast is by double passage through the 
same SLM via retroflection or phase conjugation^. In principle, if the contrast ratio for 
signal passage is C, the contrast ratio for double pasage will become C . The experi¬ 
mental configuration used to verify this fact is shown schematicaly in Fig. 23. An 
expanded and collimated laser beam from an argon laser (514.5 nm) is transmitted 
through a circular aperture (diameter = 14 mm) located at the input plane. Half of the 
aperture is masked by a NDF to form a binary amplitude mask at the input plane. The 
transmitted beam is collected by a lens (L) and directed into a barium titanate crystal. 
The crystal is oriented so that the incoming beam is phase-conjugated back via self- 
pumped mode. After reversing through the aperture, the phase-conjugated beam is 
sampled by a beam splitter through imaging optics which reimages the binary mask on 
the output plane. The input intensity distribution (i.e., single-passage through the mask) 
and the output intensity distribution (i.e., double-passage through the mask) are measured 
by scanning a detector with a small aperture (diameter = 1 mm) across the beam 
diameter at the input and the output planes, respectively. The experimental results (for 
a mask with optical density OD - 0.5) are shown in Fig. 24(a) and (b). Apart from some 
imperfection at the edge of the output image, the output contrast ratio is approximately 
the square of that of the input as expected. Instead of a phase-conjugate mirror, one can 
also place a plane mirror at the input plane immediately after the mask to retroreflect 
the beam back through the binary mask. The resulting intensity distribution is shown in 
Fig. 24(c). The phase conjugation has the key advantage of being intrinsically self- 
aligned. The grating formation time for self-pumped phase conjugation using a typical 
barium titanate crystal at a few tens of milliwatt of optical pow'er, is, however, 
relatively slow (a fraction of a second to several seconds). Note that the output obtained 
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lag. 24 Intensity profiles and contrast of a binary image: (a) at the input plane, (b) at 
tin' output plane after double passage via phase conjugate reflection, and (c) at 
’he output plane after double passage via mirror refection. 

bv ini"''» : ••!lion i , fairly uniform (see lag. 24(e)) and the rel lection is practical!* 
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I he pi k'Iih i'! ; m live iiiUti 01 mec tion holograms used to diffract the (jump beam 
iii’.a the jiia! c ban no Is spei all “d bv the matrix mask do not degrade tile contrast 
p! o . idi’J that the photoretractive crystal is placed at the Fourier plane (see Appendix 
5..S). In comparison, photoretractive image amplification via two-wave mixing in the 
image plane tends to decrease the contrast ot the input image. An experimental result 
vantv irn: this contrast preserving proper l> is given in lag. 25. Intensitv profiles of a 
hitiui v image and its amplified output are shown in lags. 25(a) and 25(b), respectively. 
Note it, it tile contrast ratio ot the input and the output image are almost : Jenin.a. and 
that !fa- vrtk'a! scale m tag. 25(h) n. 2a times that of lag. 25(a). 


SCS001 7 



(a) 


(b) 


tag. 25 I.Aperimental results showing contrast preserving image amplification by 
photorefractive two-wave mixing at the Fourier domain; (a) input intensity 
prol ile, and (b) amplil ted output intensitv profile. 


F l-.erl on the disc ussion given above, one can extrapolate that if an SLM ^.ith a 
ft: ■ - ii i o i; i he or der of 100 (tor example, a 1 error* lei trie liquid < r vatu I SLM) : used 
: a . .age mod' 1 , the macauum number ol channels (M) ot ai. M • M iros,bar 
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3.5.7 Experimental Demon strati on of Reconfigurability 

Using an experimental configuration similar to the one illustrated in Fig. 20, 
with the photographic mask (transparency) replaced by an SLM, we have demonstrated 
the reconfigurability of the interconnection by writing new interconnection matrix 
pattern on the LCTV. The pattern is generated by an IBM PC and sent to the LCTV. 
With optical power of the order of tens of milliwatts, the frame rate is currently limited 
by tr.e photorefractive response time to a few frames per second. As explained in 
Section 3.5.3, the data rate is not limited by the photorefractive response time and can 
be higher than several megahertz. To enhance the contrast of the binary pattern, the 
signal beam is passed through the SLM twice via a retroreflecting mirror (as described in 
Section 3.5.6) prior to interacting with the pump beam in a photorefractive barium 
titanate cm sta!. 
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Optical Matrix-Vector Multiplication Through Four-Wave Mixing 
in Photorefractive Media 
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Optical matrix-vector multiplication through four-wave mixing 

in photorefractive media 
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We propose and describe a new. method of optical matrix-vector multiplication by using four-wave mixing in 
photoretravlive media I'sing a BaTiO crystal, we have demonstrated such a parallel multiplication The results 
are presented and discussed. 


A large number of signal- and image-processing algo¬ 
rithms can be expressed in terms of matrix operations. 
The multiplication of two matrices is one of the most 
basic operations in matrix algebra. Optics, with its 
inherent parallelism, can offer great improvement in 
the speed of these operations. Optical processors for 
multiplying two matrices have been described in the 
literature. 1 Considerable work has also been report¬ 
ed on performing optical matrix-vector multiplica¬ 
tion. In this Letter, we describe a new method of 
performing matrix-vector multiplication by using op¬ 
tical four-wave mixing in nonlinear media. In addi¬ 
tion. we also report the first experimental demonstra¬ 
tion to our knowledge of such an optical matrix proces¬ 
sor by using a photorefractive BaTiO, crystal. 

Optical four-wave mixing has been a subject of con¬ 
siderable interest during the past several years. Much 
attention has been focused on wave-front-correction*' * 
and phase-conjugate interferometry.* 1 ” In the area of 
signal processing, optical four-wave mixing has been 
used to perform spatial information processing. 1: im¬ 
age subtraction. 1 - ■ and logic operations. 1 ' However, 
little attention has been paid to the use of the inher¬ 
ently parallel multiplication nature of four-wave mix¬ 
ing for optical matrix multiplication. In this Letter 
we describe a new method of performing matrix-vec¬ 
tor multiplication using four-wave mixing in nonlinear 
media. 

In what follows, we briefly describe the basic princi¬ 
ples of real-time matrix multiplication using optical 
four-wave mixing in nonlinear media. For the sake of 
convenience, we will limit ourselves to the case of 
square matrices, although extension to nonsquare ma¬ 
trices will be straightforward. 

Referring to Fig. 1. we consider a scheme that is 
suitable for a matrix-vector multiplication. Here, as 
an example, let us consider a discrete case in which we 
need to carry out the multiplication of an .Y-element 
vector and an .V X ,Y matrix. In this scheme, the 
vector is fanned out into .V rows of identical vectors. 
These .V x .V beamlets are directed to a nonlinear 
medium. The matrix, which also contains ,Y X ,Y 
beamlets. is also directed to the medium in such a way 
that each beamlet of tht matrix is counterpropagating 
in a direction relative to the corresponding beamlet of 


the vector. Thus, in the medium, there are .V X X 
spatially separated regions, each of which is pumped 
by a pair of counterpropagating beamlets. Now a 
probing beam, which consists of .Y X X beamlets, is 
directed into the medium in such a way that each 
probe beamlet will propagate through the correspond¬ 
ing intersection region, as shown in Fig. 1. The probe 
beam will be plane-wave beamlets propagating in par¬ 
allel. All the beamlets in this probe beam are of equal 
intensity. As a result of the four-wave mixing, each 
probe beamlet will generate a phase-conjugated beam- 
let. which, within a proportional factor, can be written 
as M(i,j)a(j). By using a cylindrical lens, a summa¬ 
tion over j can be obtained. Thus we have 

b(i) = \ Mu.j) X atj'l. (1) 

where a (j ) is the jth element of the vector a and Mii.j) 
is the matrix element. Such a scheme for matrix- 
vector multiplication can also be used for matrix-ma¬ 
trix multiplication by decomposing a matrix into col¬ 
umn vectors and then multiplying the matrix with 
each of the column vectors. 

The probe beam can also simply be a uniform plane 
wave, without beamlets. The phase-conjugated hppm 
will consist of X X X beamlets, because only a matrix 
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Fig. 2. Experimental setup used tu demonstrate the optical 
matrix-vector multiplication. 


of N X N regions is pumped by counterpropagating 
beamlets. The matrix and the vector information can 
actually be carried on any two of the three incident 
beams in the four-wave mixing process. 

The summation process can also be obtained with¬ 
out the external cylindrical lens by using a different 
scheme, described as follows. Consider a probe beam 
that consists of a column of N equal beamlets. The 
probe is incident into the nonlinear medium in such a 
way that each of the probe beamlets is made to propa¬ 
gate through a row of A’ intersection regions. The 
phase-conjugation process automatically performs the 
summation as well as the multiplication. The phase- 
conjugated beam is thus the product of the matrix- 
vector multiplication. 

As a result of the nature of the four-wave mixing 
process in the medium, this new matrix-vector multi¬ 
plier operates on the field amplitudes and thus can be 
used to handle matrices and vectors with complex 
elements. It is a coherent device rather than an inco¬ 
herent one. This aspect is distinctly different from 
most of the earlier approaches, which are all incoher¬ 
ent. When the device is operated in the coherent 
mode, the phase of each beamlet must be maintained 
uniformly over the transverse dimension of the beam- 
let. In addition, such phases must also be maintained 
Fixed in the summation process. If these phases are 
not uniform over the beamlets, the final step becomes 
an incoherent summation as a result of the spatial 
averaging. Under such circumstances, this matrix- 
vector multiplier operates on the intensities and thus 
handles only positive numbers. 

The matrix-vector multiplication is demonstrated 
experimentally by using four-wave mixing in a photo- 
refractive BaTiO - crystal. 1 ' The crystal is cut in such 
a way that the faces are all perpendicular to the princi¬ 
pal axes. All the beams are polarized extraordinary in 
the xz plane and are incident onto the a faces of the 
crystal. 

In our experiments, an Ar-ion laser beam at wave¬ 
length 514.5 nm with an output power of a few hun¬ 
dred milliwatts is used as the light source. The laser 
beam is expanded and collimated into a beam size of 
approximately 1 cm Figure 2 shows our experimental 
setup. The expanded beam is then split into three 
beams by using beam splitters BSi and BS;. To dem¬ 


onstrate the principle of operation, we chose the fol¬ 
lowing matrix and vector: 



“l 

1 

1 

0 

1 


V 


1 

1 

0 

0 

1 


0 

M = 

0 

1 

0 

0 

1 

’ a = 

1 


0 

1 

1 

0 

0 


1 


_1 

1 

1 

1 

1_ 


_i_ 


The matrix and the vector information is imprinted 
onto two of the beams by using transparencies (MA) 
and MA; in Fig. 2) that consist of circular dots. White 
circular dots represent the l’s, and the dark regions 
the 0’s. Instead of using a cylindrical lens to fan out 
the vector into five identical rows as described earlier, 
we simply use a transparency of a matrix that consists 
of five identical rows of (1 0 1 1 1). The incident 
matrices, corresponding to the matrix M and the vec¬ 
tor a, are showm in Figs. 3(a) and 3(b). Four lenses, Li, 
L_, L3, and L 4 (with focal lengths /j = f; = 0.5 m and /_ 
- fi = 0.2 m) are used to image the transparencies MAj 
and MA; onto the center of the nonlinear crystal, 
where the diameter of each beamlet is 0.1© mm with a 
center-to-center distance of 0.51 mm. The size of the 
whole image inside the crystal is 2.2 mm X 2.2 mm. 
The uniform beam is also reduced (and collimated) to 
a diameter of 3.6 mm inside the crystal by lenses L=, 
and Lr, (with focal lengths /,- = 0.6 m and /* = 0.25 m, 
respectively). The angle 9 between the uniform beam 
and one of the pump beams is 11°. The two beams 
that carry the matrix and the vector information, re¬ 
spectively, are incident onto the opposite a faces of a 
photorefractive BaTiCU crystal. The incident angle ti> 
is approximately 20° with respect to the surface nor¬ 
mal. As a result of the four-wave mixing, a phase- 
conjugated beam, consisting of a 5 X 5 beamlet pat¬ 
tern, is generated and is shown in Fig. 3(c). A cylin¬ 
drical lens is then used to perform the summation. 
The resulting product vector is shown in Fig. 3(d). By 
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Fig. 3. Photographs showing the matrix and the vectors: 
(ai input matrix; (b) input vector, expanded into a matrix 
with five identical rows; Ic) resultant phase-conjugate reflec¬ 
tion iii< trix; (d) resultant vector of a matrix-vector multipli¬ 
cation after summation of to by a cylindrical lens. 
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measuring the intensity of the individual dots, we 
found that the product vector is (3 2 0 1 4), in good 
agreement with our prediction; in other words, 
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The intensities of the product vector elements are 
measured by integrating over the cross section of each 
beamlet. The experimental error is less than 5%. 
This small error is due to environmental perturbations 
and the slow response of the photorefractive crystal. 
Although our experiment demonstrates the operation 
for matrices and vectors with binary numbers (1, 0), 
this matrix-vector multiplier can also be operated in 
the analog mode to handle elements with gray levels. 
Such an analog operation in photorefractive crystals 
requires that the probe-beam intensity be weak so that 
the error due to the intensity denominator in the index 
modulation can be neglected. This material-related 
error disappears when the matrix and the vector con¬ 
tain only binary numbers (1,0). 

In conclusion, we have proposed and demonstrated 
a new method of matrix operation by using optical 
four-wave mixing in nonlinear media. In particular, 
we have demonstrated experimentally matrix-vector 
multiplication in a photorefractive BaTiO< crystal. 

The authors acknowledge helpful discussions with 
M. D. Ewbank and M. Khoshnevisan (Rockwell Inter¬ 
national) and J. Neff (Defense Advanced Research 
Projects Agency). 
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ABSTRACT 

The concept of using color-multiplexed four-wave mixing to achieve optical matrix-matrix 
multiplication with N 3 parallelism is described. Experimental results, for the case of 
2*2 matrices, using two-color four-wave mixing in a photorefractive crystal are presented. 

1. INTRODUCTION 

In this paper we propose and demonstrate the concept of color-multiplexed four-wave 
mixing to achieve optical matrix-matrix multiplication (MMM) with N 3 parallelism. This is an 
extension of our earlier work on optical matrix-vector multiplication (MVM) 1 using rhotore- 
fractive four-wave mixing. The additional degree of parallelism required by MMM (N 3 for MMM 
versus N 2 fcr MVM) is provided by color-multplexing the nonlinear wave-mixing process. In 
Section 2 the background for optical matrix processors is briefly reviewed, and the general 
idea of using color-multiplexing to convert a MVM scheme into MMM is described. In Sec¬ 
tion 3 we present a specific approach to apply the color multiplexing technique to photore- 
fractive four-wave mixing for MMM with full parallelism. The experimental details and re¬ 
sults are discussed in Section 4. In Section 5 we summarize our results and discuss some of 
the practical difficulties and limitations, as well as some advantages of this approach. 

2 ■ BACKGROUND 

Matrix algebra has been used extensively as a convenient mathematical language to formu¬ 
late the problems frequently encountered in signal and image processing and numeric 
computing. Some of the more recent applications in optical computing include optical 
interconnections, 2 optical neural network processors, 3 and optical artificial intelli¬ 
gence. 5 In particular, two-dimensional format of a matrix can be conveniently represented 
as a spatial intensity distribution of an optical beaut over a plane. In this context, a 
matrix can be viewed as a two-dimensional image and can be transferred to any designated 
plane and scaled to appropriate size by conventional optical technique. The two 
computational primitives involved in matrix algebra, namely multiplication and addition, can 
also be carried out optically using simple reflective, refractive, or diffractive optics or 
their combinations. Various aspects of optical matrix processors have been discussed fairly 
extensively in some recent review articles. 6-9 In the following, we describe briefly only 
those directly relevant to this work. 

A basic architecture for an incoherent optical MVM scheme for performing discrete 
Fourier transforms was reported by Goodman et al. 10 A schematic diagram of this architec¬ 
ture fcr multiplying a 4 * 4 matrix and a 4 * 1 vector is illustrated in Fig. 1. In 
general, a linear array of "N" sources (light emitting diodes or semiconductor lasers, for 
instance) is used to represent the input vector (dimension = N); the value of each element 
of the input vector is encoded as the optical power of each of the source. The source array 
can be addressed electronically in parallel. A transparency or a spatial light modulator 
(SLM) is usee to encode the input matrix with the value of each matrix element represented 
by the (local) optical intensity transmittance at the corresponding location on the matrix 
mask. The ayyiopriate pixel-by-pixel multiplication as prescribed by the definition of MVM 
is achieved by projecting the source array to match the matrix mask, using anamorphic optics 
(omitted in Fig. 1 for the sake of clarity). The required summations are then carried out 
by another set of anamorphic optics (also omitted in Fig.l) that projects the transmitted 
light onto a linear detector array. The output of the detector represents the product. 

When the sources are addressed in parallel, the two dimensional matrix mask allows one to 
carry out the MVM operation with full parallelism (of order N 2 ). 

The computational complexity for the multiplication of two N * H matrices is of the 
order N 3 . In other words, it requires N 3 multiplications and N 3 summations to multiply two 
H * t; matrices. The problem can also be decomposed into N sub-problems of MVM. This is 
illustrated in Fig. 2 for the special case of N = 3. Since the columns with all the ele¬ 
ments equal to zero can be ignored, each sub-problem essentially involves only one MVM. 

Using the configuration shown in Fig. 1, one can, therefore, complete the multiplication of 
two N * M matrices Ml and M 2 in N clock cycles by clocking into the source array a sequence 
^f N vectors (i = 1 to N) which together constitute the matrix Ml; the other matrix M2 is 
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Figure 1. The incoherent natrix-vector multiplication scheme; 

the cylindrical optics are omitted for clarity. 
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Figure 2. Decomposition of a problem of matrix-matrix multiplication 
into the sub-problems of matrix-vector multiplication. 

represented by the two-dimensional mask as described earlier. In this approach, . storage 
device is required at the output end since the final result, which takes the form of H 
column vectors, is obtained sequentially in H steps. In addition to the time-multiplexing 
approach with K 3 parallelism described above, one can also use a color (wavelength) multi¬ 
plexing scheme to achieve MMM with full (N 3 ) parallelism. 

A schematic diagram illustrating the basic principle of the color multiplexing MKM 
scheme is shown in Fig. 3. Each constituent column vector of the matrix Mj is color 
encoded by illuminating it with a unique color at a unique input angle. Chromatic (as well 
as angular) multiplexing is achieved by combining the output color beam strips into a single 
column with the aid of color dispersion in a prism. Multiplication in parallel of each 
color encoded vector with matrix K? is accomplished as described earlier (see Fig. 1). 

The output, after chromatic (and angular) demultiplexing via another prism, is resolved into 
the product matrix which can be image onto a two-dimensional array of detectors. In 
practice, one may have to calibrate the output to correct for the"spectral non-uniformities 
of the sour-e, the matrix mask and the detectors. 

The color multiplexing technique is by no mean limited to the incoherent matrix process¬ 
ing scheme described above. The incoherent approach is chosen as an example in this section 
to explain the basic idea of color multiplexed matrix processors mainly because of its 
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m 3 



M3 = M^M 3 

Figure 3. A color-multiplexing scheme for incoherent 
optical matrix-matrix multiplication. 

simplicity. It also serves as a reference for comparison with the coherent processor using 
photorefractive four-wave mixing to be described in the next section. 

3. OPTICAL MATRIX PROCESSORS BASED ON NONLINEAR FOUR-WAVE MIXING 

The investigations of nonlinear optical four-wave mixing for optical phase conjugation 
and spatial convolution/correlation were reported 11 '! 2 several years ago. Recently, the 
applications of nonlinear four-wave mixing for other optical processings such as MVH, 1 
parallel logic operations, 13 optical pattern recognization for symbolic substitut ion, 1 *• and 
residue-arithmetic computation 15 have also been demonstrated. Most of these applications 
utilize the same basic multiplicative properties of four wave mixing. They differ only in 
the input spatial patterns and in their particular uses of lenses, some for imaging and 
others for spatial Fourier transform. A configuration for MVM is illustrated schematically 
in Fig. 4. A matrix-mask M is illuminated with a uniform coherent beam and imaged onto the 
nonlinear medium from one side. The mask representing the vector (ai, a2, a3) is 
illuminated from the other side. The transmitted beams are expanded by cylindrical optics, 
into three identical row vectors which align pixel-by-pixel with the matrix M. In practice 


•3 



Figure 4. A schematic diagram illustrating the basic concept 
of optical matrix-vector multiplication using 
nonlinear four-wave mixing. 
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one can avoid L*ie cylindrical optics by using a two-dimension mask to encode directly three 
(or ti, in general) identical rows of the input vector. This allows one to use identical 
imaging optics in both input ports and thus achieve better beam-overlap. Inside the 
nonlinear crystal, there are nine pairs of counter-propagating beamlets, representing the 
elements of the matrix M and the vector a. When a uniform probe beam is injected into the 
nonlinear crystal, optical phase conjugation takes place via the four-wave mixing process. 
Ondet appropriate conditions, the phase-conjugate output intensity is proportional to the 
products of the three mnut intensities. 16 A beam splitter is used to sample the phase- 
conjugate output, which consists of nine beamlets representing the pixel-by-pixel product of 
the matrix and the vector. A cylindrical lens is then used to perform the row-wise addi¬ 
tions of the beamlets to obtain the product vector (b 1 , b2, bj) at the output plane. The 
experimental detail and results demonstrating the multiplication of a binary matrix (5 * 5) 
and a binary vector using a photorefractive barium titanate crystal as the nonlinear medium 
were reported recently. 1 

Experimental observation of simultaneous optical phase conjugation of several colors of 
input from a single barium titanate crystal has been reported by Chang et al. 17 Such multi¬ 
color phase corrugation via four-wave mixing in photorefractive crystal can be employed to 
perform MMM. Figure S shows a schematic diagram of such a matrix processor. This configu¬ 
ration is basically the same as the one shown in Fig. 4 except for the addition of two dis¬ 
persive elements, one in the input for color-multiplexing, and another in the output for 
color-demultiplexing. Notice that the input matrix A is decomposed into its constituent 
vectors encoded with different colors as described in the previous section. The beam illu¬ 
minating the other matrix M and the probe beam consists of all the colors uniformly distri¬ 
buted over the beam profiles. As a result of four-wave mixing of each of the spectral 
component independently, pixel-by—pixel multiplication with N 3 parallelism is achieved. The 
signal to noise ratio of the output depends critically on the the pixel inteaction, and 
hence the beam alignment of each of the spectral components inside the crystal. In the 
scheme that we propose (see Fig. 5), the alianment is accomplished in two steps. The input 
angles for each of* the spectral components and the associated vector masks are first adjus¬ 
ted so that, the output appears as a single column consisting of all the spectral components. 
Cylindrical optics are then used to expand this single column to coincide with the matrix 
inside the nonlinear medium. This procedure ensures that all the spectral components are in 
unison and can be aligned simultaneously in the second step. Cross-talk due to interactions 
among different colors does not exist, provided the frequency separation between the colors 
are far bevonc the temporal (frequency) bandwidth of the nonlinear process. This is usually 
the case for most of the photorefractive materials since their response tine typically 
ranges from micro-second to second. 


sas ujiu 



Figure 5. A schematic diagram illustrating the basic idea of 
optical matrix-matrix multiplication using color- 
multiplexed four-wave mixing. 
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4. EXPERIMENTAL CONFIGURATION AND RESULTS 


Using a photorefractive barium titanate crystal as the nonlinear medium and an argon ion 
laser, we have demonstrated the multiplication of 2 * 2 binary matrices. The experimental 
configuration is shown schematically in Fig. 6. Broad band reflectors are used for the 
laser resonator so that it oscillates at more than one wavelength. Specifically, we operate 
the laser (Spectra Physics Model 171) at five wavelengths (514.5 nm, 496.5 nm, 486 nm, 476.5 
run, and 457.9 run) in the blue-green region of the visible spectrum. The 5-color output beam 
is expanded to about 1 cm in diameter and collimated by the beam expander/collimator (BE). 

A neutral density filter (NDF) with OD = 1 is used to reduce the phase conjugate feedback 
into the laser, and a polarizing beam splitter cube (PBS) is used to filter out the s- 
polarization component. The beam splitter BSi (R = 50%), lens Li (f = 75 cm), mirror M 2 and 
lens L 2 (f = 33.3 cm) constitute the optics for the probe beam. A 2 x 2 matrix mask MMi is 
illuminated with one of the pumping beams split off from the beam splitter BS2. The folding 
mirror M2, and the two lens imaging system llenses L3 (f = 50 cm) and Lv (f = 30 cm)) are 
used to project a reduced image of the mask MMl onto the crystal. The angle 8 subtended by 
the two beams is about 7°. The beam transmitted through BS2 is sent through a prism Pi to 
angularly separate all the 5-color components two of which (476.5 nm and 514.5 nm) are 
selected to illuminate the masks MM2 and MM3, respectively. The two color-encoded patterns 
are combined at beam splitter BS3 and imaged onto the crystal through the folding mirror M7, 
and the two lenses Ls (f = 50 cm) and L6 (f = 30 cm). For proper operation, beam alignment 
and pixel interaction inside the crystal are critical. The barium titanate crystal (5 mm x 
4 mm x 4 mm) is cut so that the angle $ between the c-axis and one of the input surface is 
about 30°. This is approximately equal to the angle between the c-axis and the grating 
(normal) vector when the incidence of the two writing beams are near normal (see Fig. 6). 
Such a configuration allows us to maximize the effective electro-optic coefficient without 
using oblique incidence, which reduces the useful area of the nonlinear medium. The phase 
conjugate beam encoded with the pixel-by-pixel multiplication of each of the colors is 
sampled by the beam splitter BS*.. In the output port, the folding mirror Ms and lens L? (f 
= 80 cm) image the center of the crystal onto the screen S, the prism P2 angularly separates 
the color components (color demultiplexing) of the output beam, and the cylindrical lens CL 
carries out the row-wise summation for each of the colors. 



Figure 6. The layout for demonstrating the multiplication of two 2x2 

matrices using color-multiplexed photore f ract l vo four-wave mixing. 

Figure 7a shows the two matrices M and A, the decomposition of matrix A into column 
vectors, and the color encoding scheme. The three masks MMi, MM 3 and MM 2 (see Fig. 6) for 
encoding the matrix M and the first and second column of the matrix A, respectively, ate 
shown in Fig. 7b. The experimental results are shown in Fig. 8. Figure 8a is a picture of 
the output pattern on the screen S when the summing cylindrical lens (CL) is renoved; the 
relative optical power of each pixel monitored by a detector is indicated by the traces from 
the chart recorder. Figure 8b is the final result representing the product matrix when the 
cylindrical lens (CL) is in place. In practice, calibration to compensate for the spectral 
characteristic of the source, masks, and detector may be necessary. 
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Figure 7. (a) The two matrices (M anc A) chosen for our demonstration, the de¬ 

composition of A into column vectors, and the encoded colors. 

(b) The three matrix masks KMi, KM 3, and MM 2 (see Fig. 6) for encoding 
• the matrix M and the two column vectors Ai and A2 derived from matrix A. 
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Figure 8. (a) the output pattern with the summing cylindrical lens (CL) (see Fig. 6) 

removed; the relative power of each pixel is indicated by the chart re¬ 
corder traces and the color of the pixel is indicated below. (b) The cor¬ 
responding results representing the product of M and A when the summing 
cylindrical lens (CL) is in place. 
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SUMMARY AND CONCLUSION 


5 . 


We have presented the basic concept and experimental results for optical matrix-matrix 
multiplication using color-multiplexed four-wave mixing in a nonlinear medium. To Clarify 
the concept of color-multiplex 1 ng, an incoherent optical approach was also introduced. 
Although the latter is relatively simple and straight forward to implement, no experimental 
results, to the best of our knowledge, have yet been reported. For multicolor four-wave 
mixing, a photorefract i ve barium titanate crystal is chosen for our experimental demonstra¬ 
tion because of its strong nonlinear response in the blue-green region of the visible spec¬ 
trum and its availability. When the number of colors is limited to a few, no significant 
reduction in phase-conjugate reflectivity of each color component is observed. However, as 
the number of colors increases, say do 10 or 20, it is anticipated that the grating con¬ 
trast, and hence the phase conjugate reflectivity for each culor will be reduced because of 
the photo-excitation of charges by all the other colors sharing the same volume. Multiple 
quantum well structures with different sets of layers designed to enhance the optical non¬ 
linearity at different wavelengths are potential candidates for color multiplexed four-wave 
mixing. With the current technology, the nonlinear optical approach is limited by the 
availability of an optical medium with sufficient nonlinearity (for four-wave mixing) over 
a wide spectral range and a multi-color coherent source that matches the same spectral 
range. 

In our experin.ental demonstration, we chose to use binary matrices and monitored the 
final result in a mixed binary format. For analog operation, several important issues such 
as the dynamic range and pixe1-to-p i xe1 uniformity require further investigation. 

One of the major advantage of the four wave mixing configuration is the availacility of 
three input ports. As a result, the system is sufficiently flexible that it car be adapted, 
with little or no modification, to perform other types of processing such as parallel logic 
operations, 13 pattern recognition for symbolic substitut ion, 1 ^ and residue arithmetic compu¬ 
tation. 15 It is quite conceivable that the incorporation of color multiplexing, vetch adds 
another dimension to t.-.is approach with even more decrees of freedom, will equip tr.e system 
for other forms of image or cate processings. 
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Optical Matrix-Matrix Multiplication with N Parallelism 
by Spatial Convolution via Four-Wave Mixing 

Arthur E. Chiou and Khanh B. Nguyen 
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(805) 373-4464 

SUMMARY 

The convolution of two (two-dimensional) spatial patterns using nonlinear opti¬ 
cal four-wave mixing in their common Fourier plane can be devised to achieve matrix- 
matrix multiplication optically, in full parallelism. In principle, it can be applied to the 
multiplication of a matrix A (of any dimension mxn) and a compatible matrix B (of 
dimension nxp). A unique feature of this approach is that the output beamlets repre¬ 
senting the elements of the product matrix lie on a straight line; therefore, a linear array 
of detectors (instead of a two-dimensional arra' ) can be used for monitoring the result. 
The major advantage of this scheme is the full parallelism achieved with a relatively 
simple optica! arrangement. The trade-off is the resulting spatial complexity. 

White and Yariv* have demonstrated that (spatial) convolution and correlation 
of two (two-dimensional) patterns can be achieved in real time by four-wave mixing in 
the common Fourier p'ane of the input patterns and recording the phase-conjugated 
output at the corresponds:, object plane. A typical experimental configuration is shown 
schematically in Fig. 1. Specifically, if U3 is a small aperture simulating a point source 
(or a delta function), a pattern representing the convolution of U1 and U2 (U1*U2) is 
observed in the output plane. An illustrative example is given in Fig. 2. The design of 
U1 and L’2 to represent two matrices so that their product is represented by U1*U2 is 
explained in the following paragraphs. Explicit mathematical analysis, including the 
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Fig. 1 A schematic diagram illustrating the basic idea of four-wave mixing in a non¬ 
linear medium located at the common Fourier plane of the input spatial 
pattern. 


set use 



U1 U2 U1 • U2 

(a) (b) (c) 


Fig. 2 A schematic illustration of convolution of two spatial patterns: (a) and (b) are 
the two input spatial patterns; (c) represents the resulting spatial convolution 
of the two inputs. 

effect of finite aperture dimensions among the others, of a similar matrix multiplication 

2 

concept based on spatial filtering was given b> R.A. Heinz et al. 






Rockwell International 

Science Center 


SC5502.FR 

For simplicity and clarity, let us consider the simplest cases of two matrices A 
and B, both of dimension 2x2, and their product C = AB given below. 


A = 

a ll 

a l 2 



a 21 

a 22 


B = 

b ll 

b l 2 



b 21 

b 22 


C = 

AB r 

c ll 

c 12 



c 21 

c 22 


(7) 


a ll b ll +a 12 b 21 a 11 b 1 2 + a l 2 b 22 

a 2i b l1 + a 22 b 21 a 21 b 21 + a 22 b 22 


A transparency U1 consists of four small apertures, each with its intensity transmittance 
proportional to each of the matrix element a- is shown in Fig. 3a. A similar trans¬ 
parency U2 corresponding to the transpose of B (i.e., rows and columns interchanged) is 
shown in Fig. 3b. Note that the vertical distance between the elements in U1 is designed 
so that it is considerably larger than that in U2, while the horizontal distances between 
the elements are identical in U1 and U2. Note also that the dashed rectangular box con¬ 
necting the four elements in U2 is an artifact to facilitate the following explanation and 
should be ignored in the actual design. By comparing Figs. 2 and 3, it is straightforward 
to see how the quadruplet of the dashed rectangular box is formed at the corresponding 
position in Fig. 3c as a result of the two-dimensional spatial convolution. As noted 
before, all the dashed lines in Fig. 3c are artifacts and the actual result should show only 
the corners of each box, i.e., a total of 16 spots with the intensity of each proportional to 
the cross product of all the elements in A and those in B. The four doublets lying along 
the Y-axis in Fig. 3c are drawn slightly offset from the axis to expose the individual 
component. In practice, the two components of each doublet are spatially overlapping on 
the Y-axis. The intensity of the four doublets, from top to bottom, on the Y-axis are 
proportional to c^i Cjj, C 22 and C 2 p respectively. The eight cross terms that do not 
contribute to the matrix-matrix multiplication can be filtered out easily as they are 
physically separated from the Y-axis. 
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C = AB 


C 11 - 8 11 b 11 + a 12 b 21 

C 12 = aii b i2 + a 12 b 22 

C 21 = *21 b 11 + ®22 b 21 



Fig. 3 A schematic diagram illustrating optical matrix-matrix multiplication by 
convolution for the case of 2 x 2 matrices: (a) matrix A and the mask U, to 
encode A; (b) matrix B and the mask li, to encode B ; (c) the spatial pattern 
resulting from convolution of U i and U^. 

Another example illustrating the case of multiplying a 2 » 3 matrix with a 
3 * 4 matrix is shown in Fig. 4. Note that all the cross terms that do not contribute to 
the matrix-matrix product are omitted in Fig. 4c for the sake of clarity. 

Using an argon ion laser (5145A) and a BaTiO-j crystal, we have experimentally 
demonstrated the concept described above for the case of 2 x 2 matrices. The images of 
the masks used to encode the two matrices A and B are shown in Figs. 5a and 5b, respec¬ 
tively. The experimental result representing the product is shown in Fig. 5c. An over¬ 
exposed version of Fig. 5c is shown in Fig. 5d to reveal the noise resulting from the other 
cross terms that do not contribute to the matrix product. 

In principle, the basic concept described above can be appied to the multiplica¬ 
tion of a matrix (of any dimension mxn) and a compatible matrix B (of dimension nxp). In 
practice, the signal-to-noise (S/N) ratio is expected to degrade as the dimension of the 
matrix increases. The problem of S/N is less serious in the mixed binary mode. The 
major advantage of this scheme is the full parallelisni ac h ieved with relativey simple 
optical arrangements. From the technological point of view, the fact that all the output 
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Fig. 4 An illustration similar to Fig. 3 for the case when the dimensions of the 
matrices A, B and C are 2 > 3, 3 * 4 and 2*4, respectively. 
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OVER EXPOSED 
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Fig. 5 Experimental results for optical matrix-matrix multiplication by convolution 
via four-wave mixing in the spatial frequency domain, (a) and (b) are the 
matrix masks used for encoding the matrices A and B ; (c) shows the result of 
their spatial convolution; (d) is an overexposed version of (c) to reveal the 
cross terms that do not contribute to the matrix product. 


elements lie on a straight line can be significant, since this allows one to use a linear 
detector array instead of a 2-D array for output detection. 
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SUMMARY 


A variety of optical matrix processing schemes [1. 2] use two-dimensional spatial 
modulation of optical intensity to represent a matrix to exploit the inherent parallel nature o: 
optics. Such an approach typically requires the projection of a spatial pattern to match another 
pattern to perform the element-by-element multiplication. The basic incoherent matrix-vector 
multiplication scheme [3], for example, requires the use of anamorphic optics to project a linear 
array of sources (or a one dimensional spatial light modulator) to precisely match a two 
dimensional matrix masks. In the matrix-vector multiplication scheme using four-wave mixing 
in nonlinear media [4], simultaneous alignment of all the pixels of the matrix and vector is a 
major task, particularly for a large number of pixels. Misalignment of the pixels may lead to 
severe errors. For a given size of matrix mask, the density of elements increases as the 
dimension (N) of the matrix increases. As a result, the requirement on alignment becomes 
more and more stringent. In practice, the critical alignment required is likely to impose a 
practical limit on the optimum dimension of the matrix (N) to be of the order of one hundred or 
less depending on the specific architecture. 

In this paper, we report a new scheme for optical matrix-vector multiplication that uses a 
phase conjugator (with a finite storage time) in conjunction with a spatial light modulator 
(SLM) to eliminate the pixel-by-pixel alignment requirement at the cost of some reduction in 
parallelism (N 2 /2 instead of N 2 ). Phase aberration due to imperfection in optics is also 
self-corrected by the phase conjugation process. The optical sj tern involved is relatively 
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simple compared with the other approaches. Without any modification, such a scheme can also 
perform matrix-matrix multiplication with N z /2 parallelism. 


LASER BEAM 



a 


Fig.l A schematic diagram illustrating the basic concept cr optical matrix-vector 
multiplication using a photorefractive phase conjugator in conjunction with a 
spatial light modulator. 


Referring to Fig.l, we use a SLM to in.press the matrix ana vector information in 
uer.ee to an input laser beam. This beam is directed toward a phase conjugator which has a 
to storage time ta photorefractive barium titan ate. for example). A cylindrical lens is inserted 


in the phase conjugate output beam path to perform the summation. 


The principle of operation is as follow, the SLM first impresses the matrix information 
onto the input laser beam. This beam is then incident into a phase conjugator which stores the 
matrix information after a finite grating formation time. When the matrix information is 
removed from the SLM, say by turning all the pixels into maximum transmission condition, the 
phase conjugate beam which contains the reconstruction of the matrix information exists for 
finite duration. 7'his finite storage time depends on the strength of the input (read) beam 
During this time, if the next frame of the SLM carries the vector information, par.-.'.he. 
mui;m;:..:tion i> performed as the phase conjugate beam propagates back through the SLM 
here tt.e vector is represented as a tssu-uimeiisjOiiu! a me. v ot N identical column vectors, v*. mm. 
v - :• -me mimi-or ot the vector. A cylindrical lens m. the output port m u .c to perform m 
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summation. The dark storage time during which the matrix information can be retrieved :> 
determined by the photorefractive material and the pumping configuration. I: ranges from 
seconds to microseconds. 

The system can also perform matrix-matrix multiplication by time multiplexing. In thus 
case, each column vector V 1 (i=l to N) which constitutes the second matrix Mi is sequentially 
impressed onto the beam to multiply with the first matrix Mj according to the matrix-vector 
multiplication scheme described above. To avoid the degradation of the information of M, 
stored in the photorefractive hologram during the readout, it is necessary to refresh the 
holographic memory with Mj to restore its diffraction efficiency. This can be done by 
re-impressing M j onto the beam after each readout cycle. Consequently, a total of 2N clock 
cycles, consisting of N cycles of write and N cycles of read, will be required to carry out the 
multiplication e: two NxN matrices. 

Using a photorefractive barium titanate crystal as a phase conjugator in conjunction with 
a -kx-S magneto-optic spatial light modulator .SIGHT-MOD SMD4SI from Semetex Conti, 
we have demonstrated the basic principle described above. Preliminary experimental results 
will be discussed. Some advantages and disadvantages of this approach will be compared with 
those of the others. 

This work is partially supported by DARPA/AFOSR contract No.F49620-S7-C-0015. 
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Optical interconnection linking laser ar'ays and detector 
arrays plays a key role in opti- al computing. 1 Conceptual¬ 
ly. such an interconnection can be achieved by using an 
optical matrix vecto, multiplication. 

v - M\. Hi 

where v is the inpu’ vector representing the signals carried by 
an arrav of .V la-er- and v' is th‘- output vector representing 
the signals received bv the array of A detectors M i> an A x 
A mairix representing th r - interconnection pattern 

When a ‘ran-pare ncx or a spot lal light modulator i SI.M 1 is 

used a- the interconnection mask, parallel matrix vector 
multiple .it n ; can be achieved I.et the interconnection 


mask consist of an .V X .V array of transmission windows. 
The input vector (a row of .V elements) is fanned out by 
using, for example, a cylindrical lens pair, so that the light 
from each element is broadcast over a corresponding column 
of the window.- After transmitting through the intercon¬ 
nection mask, similar optics are used to collect the light from 
each row of windows and to sum the'output into a column 
vector of .V elements. Such an architecture provides the .V- 
parallelism. However, a large fraction of energy is absorbed 
by the transparency or SLM. When used as a crossbar 
switch, such an architecture has an energy efficiency of only 
P'.Y. where .V is the dimension of the array. ’ This occurs 
because (X - 1) '.Y of the light energy from each element of 
the input vector does not pass through the crossbar mask. 
The energv loss increases as the dimension of the array 
increases. This is sometimes referred to as the fanout ener¬ 
gy loss. For a 1000 X 1000 crossbar switch, the loss due to 
fanout can be as high as 99.9^ This is not acceptable in high 
speed computing because signals are passing through the 
SLM at gigabit rates and the energy lnss can be enormous. 
In addition to the inherent tannin energy loss, all SLMs have 
finite insertion loss due to imperfect transmissions such as 
(corn absorption and scattering ll w e include such insert ion 
losse>. the energy efficient-y for a crossbar switch would 
bet "me t X. where t is the transmittance (f < l 1 of each 
window 
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Fig 1 Schematic drawing of an 1 x .V optica! interconnection 
u>:m: dvnamic photoretractive hologram- with A = -1 



K.g . Sc h'-mat i. draw.rig “t an \ x A npn.a’ iiiieroiiainti'-. 
u*:ng ;i«;m phot. .-etr.n :iv»- hologram- with * ; 


It h.i' been know it Cur si -me time that a h> >l> ‘graphic opt wal 
element 'HOKi can be employed tor tree space optical inter¬ 
connection. Light iront each laser source oC the input array 
is Bragg scattered and redirected to one or more detectors ni 
the output array. There are several specific requirement- 
that need to he met for an HOF. to he used for interooniu-c 
ti--n of VLSI circuit'. The-e include near perfect alignment 
and high detraction efficiency. In addition, a new hologram 
i- needed tor each new interconnection pattern. In view of 
these issues. ;• i- desirable to use dynamic holographs media 
sue h as photi.retractive crystal- ‘ Photoretractive crystai.- 
sui h a.- BaTiO. and SBX are h\ tar the most efficient real- 
time holographic meri.a. The re. rdmg ot a volume holo¬ 
gram inside these cry-t d- require-... tune ot t he order of 1 m- 
using light inten. :ty of the order of 1 W i m-.' 

In th:-. <.mmunit atui.n. we prop. ,se a st heme of forming a 
uynarnii hologram inside a phop.rrtractive crystal so that 
the hologram is capaiiie ot performing ret onfigurable optica! 
interconnection with a very high energv efficiency Such a 
sc heme utilizes the nonrecinrocal energy transfe- in two- 
wave mixing to achiev e an extremely hi: h energv efficient y 

Heterrmg to Fig- i and 2. we describe thi- new method of 
ret onfigurable optital interconnection which utilize- the 
n .nreciproca! energv tran-Ier in photo refractive two-wave 
mixing to improve the energv efficiency. Figure 1 describes 
a 1 - If ca«e for the sake of clarity in explaining the concept A 
small ! rat t ion of t hi - incident la-er beam is coupled out of the 
beam by u-mg a bean, -plitter < BS ■ Thi- sni;. 1 ! traction 
1 1 anod probe beam, i- then expanded hv using a cylindrical 
len- and pas-t - thr .ugh t spatial light modulator tSl.M.. 
Ir tin * xmnpie -fiown the la-er heap, i- to be connected to 
dele- tor- ' and a a - prc-i. ribed I.v the SI.M The t rail-la it - 
ted probe beam i- then re.. . .tiibme-i with the main bean, 
in-ide a photoretrai live i rv-tal A - a re-iilt of nniia iproca! 
ett.-rgv - oupfing aim-all the energv in the mam beam i- 
t r u. - ferret; t. - r ri* | - r. -t.i b* ill. what; came- the i r.t er. ■ - • 
net tin) patter; I j.e re-:]' i- .ir. i ;.'; . a; mien o).!,i ■ 1 1 • • 

?f 'm tP r V:j 0 P ' , C'; -. 0 Z; V, 1! 


Figure 2 describes the reconfigurable interconnection for 
la-er arrays and detector arrays. In the example (a 4 X 4 
interconnection i shown, laser 1 is to be connected to detec ¬ 
tors b and c. laser 2 is to be connected to detector- u and <1. 
laser 3 is to be connected to detectors c and d. and laser 4 is to 
be connected to detectors a and c. In terms of the matrix- 
vector multiplication, such an interconection can be written 
a- 



where c .. cc,. and i 4 are signals carried by the laser- 
element- of the array c, 

A cy lindrical len- is used to focus the 2-1) array of beam- 
into a vector 11-1) arrayi. As a result, detect..r a receives 
signals from lasers 2 and 4. detector b receive- signal- irom 
la-er 1. detector < receives signals from lasers 1and 4. and 
detector a receives signal- from laser- 2 and Such a 
concept can be extended to interconnect .V laser.- with M 
detector- where .V and M are two large numbers. 

I he two-wave mixing describee; in. Fig- 1 and 2 mav be 
viewed a- a real-time holography in which the recording and 
readout occur simultaneously inside the photoretractive 
cry.-tai. The beam splitter and the- SI.M are u-ed to rec ord a 
volume hologram which represents the interconnection p.-.t- 
tern as prescribed by the SI.M. The energv coupling in¬ 
volved in the two-wave mixing ensure- that the detraction 
efficiency during the readout i- almost ltio ,. This, ot 
course, requires a proper orientation ot the crystal so that the 
energy of the readout beam is greatly depleted. The high 
energy efficiency result- from the fact that most of the energy 
is carried by the readout beam, which does not pass through 
the SI.M but i- diffracted into the interconnect!! n pattern 
by the hologram. 

'1 he energv efficiency of such an interconnection pattern 
can be estimated tor a crossbar switch a- follow- Let ft be 
the reflectance of the beam splitter. It is legitimate to a-- 
sume that the beam splitter is practically lo—ie-.-. We may 
ai.-o assume that the surface of the photoretractive crystal is 
antireflection-coated so that the Fre-nel reflection loss can 
be neglected. In these conditions, the two beam- tha. arrive 
at the photoretractive crystal have energie-' 1 1 - fti and ft; 
.V. respectively . Inside the crystal, these two beams undergo 
photoretractive coupling. As a result, most of the energv of 
the pump beam (1 - fti is transferred tothe probe beam ft; 
A. which contains the interconnection pattern The energy 
efficiency can he easily derived and is given by 

IR 1 * m 

r - exp' — .ii. i. . ... 

A 1 + rr. exp' 

where rr, is the beam intensity ratio. 

i; - ft. v 


and I. is the interaction length, y is the coupling constant, 
and o is the bulk absorption coef fie lent. For photo retract iv e 
crystal.- such a- BaTiO and SBN. the coupling constant is 
very large u.e . •>/. >> le The efficiencv can in written 
approximate lv 


f;ft 

La 



\S e notice that tor large A . the energv efliciem v i- hmi'e-d lw 
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to amplif.ed probe intensity (vertical sale = 50mv div honzonu! scale = lgsee div 

Fit :< ( a i Signal carried In the fir..be bean, Vertical -t ali = 
mV di\. h- ruor.ta! m a <- = 1 u» d:\ > l> Amplified -izna! car re- : by 

t he pr: iite bear:, u!: er ■jn.iert 1 .’.mg tw "-w a\ e mixing Vertical -i ate = 
r»> if.V dr., hi.rizor.: %. - alt = 1 u> div 


the try-ta! bulk ;i:>' • rf>! i..n exp’ I. • anti i> maximized by 
u-mg a beam 'putter with a very small reflectance h ii.e.. /•«' - 

tn. 

Ir. car preliminary experimental work. we invest istated t he 
issue ul energy eti'iciem\ and the capability of nigh data rate 
transrni:—i< n. In examining the energy efficiency, we u*ed a 
laser beam (rum an argon-ion la*er operating at 514.A nrn. 
The 1 a'er beat:. is collimated into a beam of g-mm diameter 
by nr,ng a len- • 4 local length / = :n,. A beam splitter witf. a 
reflectance oi H - n.n.'i ir used to redirect of the energy 
through the spatial light modulator, estabi’shing the pn.be 
beam. The remainder of the energy Iran- nits through the 
beam 'phtter and const it me.' the pump beam.. In our pre 
liminarv experiment, the SLM was replaced with a neutral 
density filter wit h a \ariable optica: density to >imulate the 
lam i i - t nergy I he two beams intersect inside a fiarium 

titanate ry-ta! with an interaction length of — .*> mm The 
fi'iiiip In air. wti an mten'iiy o: p I \\ < n■■ enters the cry-tu! 
a: near m rm,,d im :dt The pr..b»- bean, with an inten»iiv 
"f o' W i rr: i' ini id’T.: a: an a:ig!> ..f 4"' I he cry-’ui i- 
or:. :.’eh - •*■,.!* the pr”bf beam amp!:;ii-d at the expeil-eof 
thi p ,mp beau. Auer p,i"i:ig through tl.i crvstal. the 
put,up a. an. I - :r: aaiic depf ’ i- i when a- the in:. :l-: t v of tin 

pr. :,e i- ep \a’. <i i.. ; i \\ , n. I ;n -at:.- *-\ penmen’ 
i • r> i'i i u-.ng in r.il <). n- w .’ r: t ran -nut Pine e ot o | 
I :.. -I I T. P. int. p-.’ • • 


0.05 and 0.005 VV/cm-. respectively. The intensity of the 
amplified probe beam remains at 1.9 \\Vcm-. This trans¬ 
lates to an energy efficiency of —20V In other words, if the 
fanout loss is 99 r c (e.g.. a 100 X 100 crossbar), the energy 
efficiency of our new scheme can be at least 20 times better 
than that of the direct approach. Notice that our crystal is 
not antireflection coated, and ~18T of energy is lost at the 
front surface. Bulk absorption in this particular crystal 
accounts for ~60 r c of the energy loss. 

fciven though the reconfiguration time is limited by the 
photorefractive response time, which is typically of the order 
of milliseconds at modest intensities, the photorefractive 
interconnection system can accept very high data rate sig¬ 
nals. To demonstrate this fact, temporal modulation was 
impressed on a laser beam (argon. ,V = 514.5 nmi using an 
acoustooptic device to simulate a signal, which is to be inter¬ 
connected with some output. The signal used was a pulse 
train of frequency f, - 0.633 MHz w ith each pulse being —0.2 
ns wide. This rate is clearly much higher than the reciprocal 
of the photorefractive response time. The modulated laser 
beam was then split into two beams and mixed in the crystal 
as described before, and the amplified probe beam was moni¬ 
tored with a photodetector. The upper oscilloscope trace in 
Fig. 3 shows the input probe signal, and the lower one shows 
tht amplified probe signal. The results show a steady-state 
response in which the temporally modulated pump and 
probe beams interact simply by di-traeting off the stationary- 
index grating that is created in the crystal alter the photort- 
fractive response time. As this experiment was performed 
merely to demonstrate the high signal bandwidth of the 
system, optimization of the parameters was not done, and 
the result- in Fig. 5 can clearly he improved. 

Summarizing: we have described a new method of recon- 
figurable optical interconnection using dynamic holograms 
in photorefractive crystals. Such a new method provides a 
very high energy efficiency. The interconnection can be 
reconfigured by using a different SLM pattern. The holo¬ 
gram formation time will limit the reconfiguration time. 
Once the hologram which contains the interconnection pat¬ 
tern is formed inside the photorefractive crystal, such a 
scheme is callable of providing the interconnection for high 
data rate transmission. In such art interconnection, the 
output of each laser can he input to any one or all the 
detectors. Optical phase conjugation can also be used in 
conjunction with the two-wave mixing to correct for any 
phase aberration that may he caused by the crystal imperfec¬ 
tions. ‘ 


The authors acknowledge helpful discussions with their 
colleagues Monte Khoshnevisan and Derek Cheung. This 
research is supported in part hv DARLA AFOSR under 
contract F4962O-87-C-O015. 
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Abstract. This paper describes various nonlinear optical phenomena in 
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cation and subtraction, logic and matrix operations, and optical intercon¬ 
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1. INTRODUCTION 

The held of optic*, speubcall) nonlinear optics, is emerging as 
an area ot increasing importance tor application* in optical com¬ 
puting because of its inherent parallelism and its extreme hand 
xxidth of communication 1 In addition, the nonlinear optical (NT.O’ 
response ot materials allow* lor the po*sibiht\ ot manipulating 


In •lieu f.i:* : ( ,f lii r.-.. c ! O.; ' !u.** re.:*eJ iii.inu*.np ; re.cixL-u Id?. 
13 ]‘v*.| tl ■ 1 t-u pi; r '!i. iln J e b po . 

1 I u*‘i S .i - * f'rv-t. - (ip-:, .i: in-injineni.it !■ >! tnetne-er- 


the propagation of an optical beam bx another beam of light 
inside a nonlinear medium. These manipulations include changes 
ot phase, polarization. intensitx. direction of propagation, and 
even the frequency of light itself."' In the area of optical com¬ 
puting. these NLO manipulations older unique capabilities in 
holographic interconnection, parallel image subtraction, phase 
conjugation, optical data storage, spatial light modulations, and 
bistable optical devices. 

Some of the most useful NLO phenomena include optical 
phase conjugation." four-waxc mixing.' and txxo-xxaxe mixing h 
Traditional!). these phenomena are observed in nonlinear media 
ib a exhibit a large third-order susceptibility. Typical!). these 
phenomena require very high optical intensities for efficient op¬ 
eration. Using one of the most efficient NLO materials such as 
carbon disulfide, an optical intensity level on the order of 1 MW 
enr is still needed These intensity levels are too high for optical 
computing applications Photorefractive materials such as BaTiO.. 
SBN. BSO. and BGO are by far the mosi efficient media for 
these NLO phenomena using relatively low intensitx levels te g . 
I W cm'i. In addition, these materials are acentric and exhibit 
nonreciprocal energy transfer in a degenerate two-beam coupling 
configuration. Such optical nonreciprocity is not available in a 
third-order nonlinearitx 

In this paper, we review some of the most useful NLO phe¬ 
nomena in photorefractive media We then discuss their appli¬ 
cations in optical computing 

2. PHOTOREFRACTIVE MATERIALS 

The photorelractive eflecl is a phenomenon in which the local 
index ot retraction is changed by a spatul variation o! (he light 
intensitx Such an effect was first reported in l ( thh The spatial 
index variation leads io distortion ot the wavefront, and such an 
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effect was first referred to as "optica! damage." The photore¬ 
fractive effect has since been observed in many electro-optic 
crystals including LiNbCh. BaTiOv. SBN, BSO. BGO. GaAs. 
InP. etc. It is generally believed that the photorefractive effect 
arises from optically generated charge earners that migrate w hen 
the crystal is exposed to a spatially varying pattern of illumi¬ 
nation with photons having sufficient energy. Migration of the 
charge earners due to drift or diffusion produces a space-charge 
separation, which then gives rise to a strong space-charge field 
Such a field induces a refractive index change via the Pockels 
effect. This simple picture of the photorefractive effect explains 
several interesting steady-state optical phenomena in these 
media. 

Althoueh there are several models for the photorefractive 
effect. K_1 ~ the Kuk; tarev Vinetskii model is the most widely 
accepted one 8> In this model, the photorefractive materials are 
assumed to contain donor and acceptor traps that anse from 
imperfections in the cry stal These traps create intermediate elec¬ 
tronic energy states in the bandgap of the insulators. When 
photons with sufficient energy are present, electronic transitions 
due to photoexcitation take place. As a result of the transitions, 
charge carriers are excited into the conduction band and the 
ionized donors become empty trap sites. The rate of carrier 
generation is (si - (BiiN'd - S'C ). whereas the rate of trap 
capture is ‘YrNNu ■ Here, s is the cross section ot photoion¬ 
ization. I is the light intensity . p is the rate ot thermal generation. 
Vr is the carrier-ionized trap recombination rate, and N and 
Nlj stand for the concentration of the carriers and ionized traps. 
No is the number density of the donor. 

Th^ ..pucc-GioJgi. field p.uducsd by tne migration of the charge 
carriers is determined by the following set ot equations’" '. 

r* r 1 

-N = -NY, - -V J . Hi 

r»t fit e 

-NV - (si - B>'N- - S-", ■ - MAN." . 

rU 

J = e) aN j F. - -y-VlocN ) - pie . Gi 

VnF i = eiNT - N - Nf. > . < J| 

where c is the unit vector along the c-axis of the crystal. N A is 
the acceptor concentration, jj. is the mobility. T is the absolute 
temperature, k is Boltzmann's constant, n is the index of re¬ 
traction. t is the dr' Irk tensor, pi is the photovoltaic current, 
e is the electronic charge. and p is the photovoltaic constant 
E stands tor the space-charge field F is the total field, which 
includes E" and any external or internal fields isuch as chemical 
or internal tem>elevtrte tieldsi 

As a result of the presence <>/ the space-charge held, a change 
in the index of retraction is induced by means of the linear 
clectro-optis effect' I Pockels el test i 

Jk(' ! I - t.K. . 

r 

where r,,i o the electro-.,pile soethcient (with l.t.k x.y.zi 

3. NON FINE. \R OPTIC AI PHENOMEN A 

Ail o: .re NI.O phenomena mentioned earlier involve the tor 
ma’ii in i-I volume in, le> crating- nouie tne ph. torelractive me 


dium. The simplest way to produce a volume index grating is 
to employ two-beam interference inside a nonlinear medium. 
This process takes advantage of the nonlinear response of the 
material to illumination bv electromagnetic radiation. Thus, when 
two beams of coherent light intersect inside a photorefractive 
medium, they create a volume index grating When the two 
beams propagate through the self-induced grating, they undergo 
Bragg scattering. One beam scatters into the other and vice versa 
in a process termed two-wave mixing. The hologram formed by 
the two-beam interference inside the photorefractive medium 
can also be erased by uniform light illumination. Thus, dynamic 
holography is possible by using photorefractive materials 

The Bragg scattering involved in two-wave mixing is very 
similar to the readout process in holography. In four-wave mix¬ 
ing. a third beam is used to read the hologram formed by the 
two-beam interference. The fourth beam is generated as a result 
of the Bragg scattering. To satisfy the Bragg condition, the third 
beam must be counterpropagating relative to one of the two 
beams that arc involved in the formation of the index grating 
In two-wave mixing, the Bragg condition is automatically sat¬ 
isfied. In w hat follow s, w e discuss some of the NLO phenomena 
in photorefractive media that are useful for optical computing 
applications. 

3.1. Two-wave mixing 

Two-wave mixing is important in practically all of the useful 
NLO phenomena observ ed in photorefractive media The process 
of torming an index v ariation pattern inside a nonlinear medium 
using two-beam interference is very similar to that of hologram 
formation. Such an index variation is often periodic and is called 
a volume grating In addition to the holographic process known 
traditionally, beam coupling occurs simultaneously. This is a 
unique property of photorefractive materials. In what follows, 
we briefly rev iew the coupling and energy transfer of two beams 
inside a photorefractive medium. 

Consider the interaction of two laser beams inside a photo- 
refractive medium (Fig 1). If the two beams are of the same 
frequency, a stationary interference pattern is formed. Let the 
electric field of the two waves be written 


t: -- Ae.xpliuut - k -rij . j = 1.2 . ifo 

where A|.A; arc the wave amplitudes, w is the angular fre¬ 
quency. and k].k; are the wave vectors For simplicity, we 
assume that both beams are polarized perpendicular to the plane 
of incidence tie. s-polarizedi 

Within a scaling factor, the intensitv can be written 

1 - L = I - I. <~i 

Using F.q th) tor the eLvtrm held, this can be written 

I - A ' - A; ' - A'A.exp: ik r - A ; A*exp'iK r> . (Si 

w here 

k k k im 

and k 2 tt \. where \ o the period ot the fringe pattern 
The intensitv [F.q (Si| has ;i spatial variation inside the photo 
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refractive medium According to Kukhtarev's model, such an 
intensity pattern will generate and redistribute photocarriers, cre¬ 
ating a space-charge held in the medium This held induces a 
volume index grating due to the Pockels effect, in general, the 
index grating will have a spatial phase shift relative to the in¬ 
terference pattern. s The index of refraction, including the fun¬ 
damental component of the intensity-induced gratings, can be 
written 


a 2 
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no 


ni . A*A; 

—exp(ub i —j—exp( 


iK-ri * c.c. 


Fig 1. Schematic diagram illustrating the basic idea of photorefrac- 
tive two-wave mixing. 


where 

— Ii -*■ 1; = Af' fA-" . (ill 

ru, is the index of refraction when no light is present. <±> is real, 
and m is a real and positive number For the sake of simplicity . 
we assume a scalar grating. The phase 6 indicates the degree 
to which the index grating is shifted spatially relative to the 
intensity pattern. In photorefractive media that operate by dif¬ 
fusion only ii.e.. no external static heldi. for example. BaTiO-,. 
the magnitude of 6 is 2 with its sign depending on the direction 
of the e-axis. K is the grating wave vector given by Eq. (9) 
The parameter m depends on the grating spacing and direction 
and on the material properties of the crystal, e g . the electro- 
optic coefficient. Expressions for niexpiid) can be found in 
Refs. 9 through 11 

The spatial phase shift between the interference pattern and 
the induced volume index grating has been known for some 
time * 1 ' The phase shift allows for nonreciprocal steady-state 
transfer of energy between the beams." To investigate the 
coupling, we substitute Eq (10) lor the index of refraction and 
E = E) * E: for the electric held into the wave equation. This 
leads to a set of coupled equation'. 

The solutions lor the intensities 1 1 1 z) and Iz) are 1 


1 - m exp'v 
I - rr. 

I-i/i = I'i();-expi-uzi . (13 i 

1 * mexp‘ - y /> 

where i is normal to the crystal surfaces and m i' the input 
mtensits ratio 


(n the absence of absorption to = Oi. 1; /1 is an increasing 
function of / and lpz) is a decreasing function r»t /. provided 
7 is positive The sign of y depends on the direction of the 
c-axis A' the result of the coupling for y > 0 in Fig 1. beam 
2 gains energy from beam I It this two-wave mixing gain is 
large enough to overcome the absorption loss, then beam 2 is 
amplified Such an amplification is responsible for the fanning, 
stimulated scattering and oscillation of iasc beams in photorc- 
fraettve crystals lh 

We mentioned the holographic interpretation of two-waxe 
mixing in photorefractive media earlier Elaborating, we con¬ 
sider the formation of art index grating due to the presence of 
two coherent laser beams inside a photorefractive crystal This 


is formally analogous to the recording process in conventional 
holography. Consider Fig. 1. in which two laser beams intersect 
and form an induced index grating. The index grating [Eq. (10l] 
contains the product of the amplitudes Ai and A:. This grating 
is a hologram formed by an "object" beam A] and a "reference" 
beam A;. The diffraction component of the transmission function 
of such a hologram [t — e.xpt - iAnL Af] is given approximately 
by 

i - An - ATA.'Cxp' — iK-r) - AiA'exptiK ri . (15i 

where A) and A; denote the complex amplitudes of the object 
and reference fields, respectively. Equation (15) assumes that 
the modulation is weak, so the higher order terms can be ne¬ 
glected. 

For reconstruction (see Fig. 1). the hologram is illuminated 
b; the reference beam A;exp' - tkyri The ditfracted beam can 
be written 

t|A ! A* A;expi - ik; r > . ilKi 

where p is the diffraction efficiency. Notice that the phase of 
A; cancels out and the diffracted beam is a reconstruction of 
the object beam Aiexpt - ikon. Similarly, the reference beam 
A; can be reconstructed by illuminating the hologram with object 
beam A|. provided beam A] is a phase object (t.e.. A| has only 
phase variation, with I.Ag = constant). 

In addiuon to holographic analogy. two-wave mixing exhibits 
amplification, which is a unique feature not available in con¬ 
ventional holography. Using these two properties, two-wave 
mixing can be used for beam processing. As a result of its real¬ 
time holographic nature, photorefractive two-wave mixing ex¬ 
hibits nonreciprocal enerev transfer without anv phase cross 
talk. 1 " 

The lack of phase cross talk can be understood in terms of 
the dtl fraction from the self-induced index grating in the pho- 
torefractive crystal. Normally, it a beam that contains phase 
information tlnr.ii is diffracted from a fixed grating, the same 
phase information appears in the diffracted beam For a selt- 
tnduced index grating, the phase information vlnr.ti is impressed 
onto the grating in such a way lhal diffraction will be accom¬ 
panied by a phase shift -ilitr.d. Such a dynamic hologram 
makes sell-cancellation of phase mtormation possible when the 
incident beam is diffracted from the grating produced hv the 
incident and reference beams This self-cancellation of phase 
tnio'niauon is actually equivalent to the reconstruction of the 
reference beam when the hologram is read out by the object 
beam 

Energy transfer without phase cross talk can be employed to 
compress both the sp.itial and temporal spectra of a light beam C 1 ' 
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Fig 2. Schematic diagram illustrating the basic idea of photorefrac- 
tive four-wave mixing 

In other words, the energy transfer without phase cross talk can 
be utilized to clean up the wavefront of a laser beam. This has 
been demonstrated experimentally using a photorefractive SBN 
crystal 19 " 1 ' to correct for the spatial-wavefront and temporal- 
wavefront aberrations. 

3.2. Phase conjugation and real-time holography 

Referring to Fig 2. consider optical four-wave mixing in a 
photorefractive medium Beam' A; and A; incident from the 
left form a volume index grating inside the medium If the third 
beam Ai is exactly counterpropagating relative to A;, then a 
fourth beam is generated thai propagates backward relative to 
A) and is a time-reversed replica of that beam. This is what 
happens in the readout process of conventional holography." 
except that in four-wave mixing the writing and readout pro¬ 
cesses occur simultaneously For this reason, phase conjugation 
bv means (if optical four-wave mixing is often referred to as 
real-time holography 

It we consider the combination of the photorefractive medium 
and the two counterpropagating beams i A; and A.) as an optica! 
device, then such a device behaves like a retroreflective mirror 
An incident beam will be phase conjugated and propagate back¬ 
ward. retracing its path Su an optica! device is called a phase- 
coniugate mirror iPC.Mi Ketiection from a PCM exhibits many 
interesting properties, including wavefront aberration correc¬ 
tion. polarization restoration. and amplification ' The am¬ 
plification is possible because energy can be coupled from the 
two counterpropagating beams. 

A class of phase conjugators that have received considerable 
attention, recently are the self-pumped phase conjugators."'* 

In these conjugators. there are no externally supplied counter¬ 
propagating pump beams Thus, no critical alignment is re¬ 
quired The phase-coniugate retlecnv ity is relatively high at even 
low laser power These devices are by far the most convenient 
PCMs available Although several models have been developed 
for self-pumped phase conjugators. - " the phenomena can be 
easilv understood with the resonator mode! In the tuo- 

interaction region model." two separate four-wave mixings are 
responsible lor the phase conjugation In the 2 K grating model, 
sell-pumped phase coniugati-m is considered as a process similar 
to stimulate ‘ Brilloum scattering 'SHS> In the resonator model, 
the crvstal k ube is viewed as an optica! cav it v tha! supports a 
multitude of modes These modes are trapped inside the crystal 
owin'-’ to total internal reflection at the surface 

When a laser beam enters the crvstal. some of the modes mav 
be excited a- a result of the strong parametric gain due to two 
wave mixing In particular, ring o-cillation can be generated 
accor-ime !• a theorv developed in Ret - ; W hen the config¬ 
uration of tne resonatice eav i’\ relative to the incident laser beam 


supports bidirectional oscillation, a phase-conjugate beam is gen¬ 
erated by means of four-wave mixing 

According to this theory. 33 the frequency of oscillation inside 
the cry stal can be slightly detuned from that of the pump beam 
If to is the frequency of the incident laser beam, the frequency 
of the internal oscillation can be written 

OJ ' = oj + b . (1 1) 

w here 8 is the frequency detuning and is of the order of ± 1 Hz 
for BaTiO?. This frequency detuning depends upon the path 
length of the nng oscillation inside the crystal The bidirectional 
oscillation provides the counterpropagating pumping beams. By 
conservation of energy, the phase-conjugated beam has a fre¬ 
quency of co * 25. 

The resonator model presents a simple explanation of the 
frequency shift observed in BaTiO? self-pumped phase conju- 
gators ' Experimental evidence indicates that internal os¬ 
cillations play a key role in the generation of phase-conjugate 
waves ' 

3.3. Phase-conjugate interferometry 

A phase-conjugate interferometer iv simply an interferometer 
that uses one or more PCM'. Several configuration' of phase- 
conjugate interferometers have been investigated These include 
a Mach-Zehnder interf erometer incorporating a pha'C conjugates 
for sell-referencing phase contouring based upon the interference 
between an incident wave and its phase-conjugate counter¬ 
pan. ' s a Michelson interferometer with one of the mirrors 
replaced by a phase conjugator. 39--1 and a Michelson interfer¬ 
ometer with both mirrors replaced by one or two phase conju- 
gators.'*'"* 9 Several optical computing applications, including 
coherent jmage subtraction. 44 ** "exclusive OR" logic opera¬ 
tion. 4- 4 and spatial differentiation. 4 have been demonstrated 
using the Michelson coniiguration with both mirrors replaced 
by a single self-pumped BaTiO? phase eonjugator (see See. 3.2 1 . 
A similar configuration for temporal differentiation tor novelty 
tilteringi has also been conceived. 49 In this section, we briefly 
review the unique features of a phase-conjugate Michelson in¬ 
terferometer (PCMIi. Applications in optical computing arc dis¬ 
cussed in Sec 4.2. 

One of the most important and interesting phenomena asso¬ 
ciated wuh PCMIs is time reversal. Referring to Fig 3<ai. con¬ 
sider a Michelson interferometer with PCMs An incident laser 
beam twith intensity 1 i is split by beamsplitter BS (assuming 
lossless) into two components, which are retroreflecied back bv 
the PCMs. According to Stokes' principle of time reversibility. " 
when the two phase-coniugate beams recombine at the beam¬ 
splitter. a time-reversed replica ot the incident beam is generated 
This beam propagates backward, retracing the original incident 
beam path. Thus, there is no light at the output port (port A), 
the total darkness is a result of time reversal Such an optical 
time reversal has been observed experimentally using two cou¬ 
pled BaTTO-, crystals for the PCMs 

Referring to Fig. 3ib). (ignoring transparencies Tqx.yi and 
T; i x. v). the efteets ot which will be discussed in See 4.2]. let 
t.r and i . r be the amplitude transmission and reflection ok’I- 
tiucnts ot beamsplitter BS tor beam entering from the right and 
let:, respectively T..c input beam with complex tick! amplitude 
I is split bv the beanisplittei mm a reflected component rl and 
a transmitted component t! T hese too components are phase 
von mealed to give r'l ' and t*|-rf and are : ■ ■ niier spin and re 
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Fig 3 (a) Time reversal for a Michelson interferometer equipped with 
PCMs (b) Quantitative analysis of the time reversal of a PCMI. 

combined with r*E*r ~ t*E’t' at port B and r*E*t * rE'r 
at port A The principle ol reversibility requires that 


(r' - t*t">E* = E* . i Ihi 

Hr* - fr'lE* = 0 . i I9i 

Thus, the amplitude coefficients satisty the relations 
r : 'i*i = I . Ooi 

tr* * t*r - 0 . !21 i 


Equation (21 1 establishes the condition tor complete cancellation 
of the held at port A. These equations are known as Stokes' 
equations 

In practice, the cancellation at port A ma> not be complete 
because one cannot obtain an ideal phase conjugation. The mea¬ 
surement of time-reversal fidelity using a PCMI was reported 
recently by hwbank and \ a/quez 1 The phase instability arising 
from the use of two independent (not coupled) self-pumped phase 
conjuealors is discussed in Ret 43 

4. OPTICAL COMPETING APPLICATIONS 

The NLO phenomena, including two- and tour-wave mixing, 
optical phase comueation. and phase coniueate interferometry. 
described in the previous section can be employed in a wide 
rune- of applications tor optical computing and image process 
me A lew examples are given in this section i" illustrate some 
ot the men:- and limitations of these NI.O approaches 


Fig 4 Experimental configuration and results of image amplifica¬ 
tion using photorefractive two-wave mixing 

4.1. Image amplification 

The energy transfer in photorefractive two-wave mixing dis¬ 
cussed in Sec. 3.1 can be used for image amplification. ' Spe¬ 
cifically. the interaction of a strong pump beam w ith a weak 
image-bearing beam inside a photorefractive medium results in 
an amplification of the image information. To preserve the con¬ 
trast of a gray scale transparency , it is crucial that the amplifi¬ 
cation factor be umlonn spatially and throughout the intensity 

dvnamic ranee of the lmaee-carrv me beam This issue has been 
* ** *• < ; *“ 
analyzed by Fainman and Lee.' Their results indicate that a 

uniform intensity gain of more than 4(1 dB over an input dy namic 

range of four decades li.e.. 13 bits) is possible in photorefractive 

crystals such as BaTiO. and SBN with strong coupling (-yL - 

If), where y is the coupling coefficient and L is the interaction 

length). 

The basic idea described above is illustrated in Fig. 4 along 
with the experimental results Both the pump and the image¬ 
carrying beams were split off from the same laser output 
(514 5 nm). The signal beam was first attenuated to be about 
six orders of magnitude weaker than the pump and then sent 
through the transparency u> carry the image information prior 
to mixing with the pump beam inside die crystal. An xyz-cut 
BaTiO: crystal was used as the nonlinear medium. The close 
resemblance between the amplified image (with a gain of ap¬ 
proximately 2().(KX)i and the original image indicate*- that the 
gain is sufficiently uniform over the dynamic range of a typical 
gray scale picture to preserve the contrast. 

In photnretracnve crvstals such as BaTiO-. and SBN. the 
phase shitt between the interlerence fringes and the index grat¬ 
ing. in the absence ot an externally applied electric held, is 
approximately rr 2 (see Sec 3 1 i As a result, the coupling con¬ 
stant and hence the signal gain arc strongest when the optical 
frequencies ot the pump and probe beams are identical |see. for 
example. F.q ( 14 1 in Ret 54; I he signal gain decreases 
drastically 1 ' as the optical frequency diticrence between the 
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pump and probe beams exceeds the reciprocal of the photore- 
fractive response time Since the typical response time of a high 
gain photorefractive crystal (e g.. BaTiOj or SBN) is of the 
order of a fraction of a second at an intensity level of the order 
of 1 W/crrT, the frequency difference of the two beams has to 
be less than 10 Hz Practically , this requires that the two beams 
originate from the same laser To overcome this limitation, there 
have been numerous investigations worldwide to search for al¬ 
ternative materials or processes w ith much faster time response. 

Several other image processing applications employ photo¬ 
refractive two-wave mixine These include incoherent to co¬ 
herent image conversion.' positive and negative image trans¬ 
fer. 50 dynamic range compression. 5 5h NLO range imaging. 50 
visualization of vibrational mode patterns. ‘ and novelty fil¬ 
ters.' 59 In Sec. 4.2 we describe a unique application in image 
subtraction using phase-conjugate interferometry 

4.2. Image subtraction and “exclusive OR" logic 
operation using phase-conjugate interferometry 

The PCMI examined in Sec. 3.3 can be readily adapted to per¬ 
form several mathematical operations over a (or a pair of) two- 
dimensional spatial patterm s) of amplitude or phase modulation 
These include parallel image subtraction and addition, negation 
and XOR logic operations, and spatial and temporal differen¬ 
tiations We briefly review some of these applications in this 
section The major advantages of the techniques using PC.Mls 
over conventional approaches are dynamic stability and align¬ 
ment insensitivity. Details are discussed in Ref. 60. 

Optical imam subtraction Optical image subtraction 0 *' ha-, 
been a subject of considerable interest because of its potential 
in many important applications. In this technique, one image is 
subtracted optically from another to detect the differences be¬ 
tween them 

The PCMI described earlier can be used as a real-time image 
subtractor Referring to Fig 3(b). consider the insertion ol a 
pair of transparencies (with intensity transmittance functions 
Trfx.y i and Trfx.v)} in the beam path in each arm of the inter¬ 
ferometer equidistant from beamsplitter BS. By an analysis sim¬ 
ilar to that in the previous section, it can be shown that when 
the two phase-conjugated beams recombine at the beamsplitter, 
the image intensity at output plane A is given by 

Uix.y > = Fi" p ' t*r T;i x.y i - r'tTax.v i • . 1 22 * 

where p is the reflection coefficient of the phase eomugators 
(assuming equal phase-conjugate reflectivity i Using Stokes 
equation [Up (21)(. this become.. 

U'x.yi = E ' p RT Tp x.y i - T;<\.v>‘ . i23i 

w here R = r" and T =• 't" are the intensity reflectance and 
transmittance, respectively, for beamsplitter B.S Note that the 
output intensity is proportional to the square of the difference 
between the inten.ity transmittance functions 

Real-time parallel image subtraction using a PC’Ml that in¬ 
corporates a self-pumped BaTiO. crystal as the PCM lat 514.5 nmi 
has been demonstrated ~ An example of the experimental 
result' is given in Fig 5 The horizontal and vertical bars are 
the linages at output port A (see f ig 3ib)j of transparencies I 
and 2 respectively, when the illuminating beam for trie other 
arm e Noshed The right checkerboard pattern represent' the 
coherent subtraction ot the two images due to destructive intet- 





Fig 5. Experimental results for image subtraction using a PCMI. The 
horizontal and vertical bars are the images of transparencies 1 and 
2, respectively, when the illuminating beam for the other arm is 
blocked The checkerboard pattern is the intensity distribution of 
the coherent subtraction of the two images. 


ference when both illuminating beams are present. Note that the 
intensity distribution where subtraction takes place is fairly um- 
torm and is very close tv' that of the true dark background ii.e.. 
dark squares where the dark regions of the bars overlap) 

"Exclusive OR" Ionic operation The image subtractor can 
aKo perform logic operation. Consider the case in which both 
transparencies arc either 1 or 0. According to Eq (23). a com¬ 
plete cancellation would require that the two transparencies be 
identical. An output intensity of 1 will appear at port A when 
only one of the transparencies transmits Thus, such an image 
subtractor can act as an XOR gate. When the transparencies are 
encoded with matrices of binary data, such an image subtractor 
acts as a two-dimensional array of XOR gates. 

Intensity inversum A special case of image subtraction is 
intensity inversion, which is obtained by removing one of the 
transparencies such that the transmittance becomes unity in one 
arm Experimental results are shown in Fig. 6 

Amplitude inutiic subtraction: The image subtraction methods 
described thus far result in subtraction of two intensity patterns 
due to double passes of light through each transparency A new 
method that provides parallel and real-time amplitude subtraction 
of two images by using holographic interference in photore¬ 
fractive media has recently been demonstrated This method 
uses a double Mach-Zehnder interferometer with a BaTiO. PCM 
and requires only a single pass of light through the transpar¬ 
encies 

Spatialdifferentiation A PCMI w ith a pair of identical spatial 
patterns (one in each arnu can be used to perform spatial dif¬ 
ferentiation Referring to Fig 3ibi. if one ot the transparencies 
is displaced relative to the other along any direction v perpen¬ 
dicular to the optisal axis by an amount Ax. Icq (23i becomes 

lev.vi - I '(i RTU.iv - Av.vi - T,.i\.yi 1 24 1 
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Fig. 6. Experimental results for the intensity inversion, (a) Intensity 
distribution of the phase-conjugate beam in the first arm with the 
transparency removed [I, *T,|x,yl 2 = 1], (b) image of the trans¬ 
parency in the second arm [l 2 1 T 2 (x,y) 2 ], and (c) the intensity in¬ 
version of (b) [I * Ti(x,y) - T 2 (x,y) 2 = 1 - T 2 (x,y) 2 ] 


where To is the transparent') intensitv transmittance This equa¬ 
tion can be interpreted as the discrete version of spatial differ¬ 
entiation with respect to the \ coordinate. Note that such an 
operation leads to edge enhancement Using a slight!) different 
experimental configuration, kwong et al. 4 demonstrated both 
first- and second-order spatial differentiation of a two-dimen¬ 
sional transmission pattern 

Temporal differentiation: An architecture that uses a two- 
dimensional spatial phase modulator in one arm of a PCM1 tor 
the temporal differentiation of a scene has been proposed" ’ and 
demonstrated/' 4 

4.3. Matrix operations 

In this section, we review the use of four-wave mixing in non¬ 
linear media to perform matrix-vector and matrix-matrix mul¬ 
tiplications. It was pointed out in Sec 3.3 that tour-wave mixing 
is an NLO process in which three input waves mix to yield a 
fourth wave In phase-matched four-wave mixing, the three input 
waves consist of two counterpropagating pump waves E; and 
£: and an arbitrary probe wave E; These three waves are cou¬ 
pled through the third-order susceptibility \" of the medium 
A fourth wave Ea is generated that can be written as' 

L, r \ 'l. ET’ 12?' 

This indicates that four-wave mixing can be used for the mul¬ 
tiplication of signals In addition, it we use the parallel nature 
of optical waves, each wave can earn spatial information for 
the purpose of image processing/’ logic operations. w '’ nu¬ 
meric processing.*' and matrix operations/" ~ In what follows, 
we describe some unique concepts that use the two transverse- 
dimensions to earn matrix information for the purpose of matrix 
multiplication 

Matrix multiplication between two N > N matrices can be 
written as 

( = AB 12n; 


w here 

( ^A,b. CT> 
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Fig 7 Schematic diagram illustrating the basic idea of optical ma¬ 
trix-matrix multiplication using nonlinear four-wave mixing 


spectivelv. Beam 1 is chosen to propagate along the y-axis and 
beam 2 along the x-axis. The matrices can he either continuous 
or discrete In the discrete case, each beam consists of a matrix 
of bcamlets. as shown in Fig 7(b) 

In the nonlinear medium, the two matrix-carrying beams form 
an interference pattern, and a volume grating is formed. This 
grating contains inlormation about the product ot the matrix 
elements of the two matrices and can be written as 

An = n.'Aix.ziB'u.v lexpdK-r ) -*- c.c . (2St 

where K is the difference ot the wave vectors of the matrix- 
carrying beams The parameter n; is the Kerr coefficient and is 
proportional to the third-order susceptibility y' " of the medium. 
Note that the nonlinear response o! the medium performs the 
function of parallel multiplication 

The volume grating is read out bv a third beam, which can 
simplv be a plane wave The diffracted beam consists ot the 
integrated contribution from each part of the grating along the 
beam path and can be written 


Vac that matrix multiplication consists ot two main operations 
a parallel multiple ate >n and a summation 

Reternnc to Fig "mi. consider a lour-wave mixing config 
uratior; suitable tor matrix multiplication Beams I and 2 contain 
(he information about the two matrices A x./- and Bi/.vi. re 


Cix.vi* J A- x./iB't/.v id.’ . ( 2^1 

where the integration is carried out along the beam path Note 
ilia! the integration completes the operation id matrix multipli¬ 
cation The information about the product ot these two matrices 
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Fig 8 Schematic diagram illustrating the basic idea of optical matrix-vector multiplication using nonlinear four-wave mixing 


is now impressed on the tranks erse spatial distribution of the 
diffracted beam. In t -actice. the dimensions of the matrices mas 
be limned by cross talk between various channels due to dif¬ 
fraction and or imperfection in the optics 

Owing to the phase-matching requirement, the readout beam 
must be incident along a direction that satisfies the Bragg dif¬ 
fraction condition to achieve high efficients In anisotropic non¬ 
linear media, the polarization states, as well as the direction of 
propagation, can be chosen such that the largest of the nonlinear 
susceptibilities is lulls utilized 

Four-wave mixing can be either degenerate or nondegenerate. 
In the degenerate case, all of the beams have the same frequency 
In the nondegenerate case, the frequencies of the beams can be 
slightly different This mas be useful for the purpose of sepa¬ 
rating the diffracted beam from the undiflracted portion. 

Another scheme, shown in Fig. 8. is suitable for matrix- 
vector multiplication Here, as an example, consider a discrete 
case in which we need to cam out the multiplication of an N- 
element vector and an N / N matrix. The vector is fanned out 
into N rows of identical vectors These N * N beamlets arc- 
directed to a nonlinear medium. The matrix, which also contains 
N v N beamlets. is also directed to the medium in such a wav 
that each beam of the matrix is counterpropagating relative to 
the corresponding beam of the vector Thus, in the medium there 
are N ' N spatially separated regions, pumped by a pairot coun- 
terpropagating beams Now. N'/N probe beams are directed 
into the medium in such a wav that each probe beam can prop¬ 
agate through an intersection region The probe beams will be 
"plane wave" beamlets propagating in parallel As a result of 
the tour wave mivng. each probe beam generates a phase-con 
lugated beam that. a ithin a proportionalitv factor, can be wruten 
M'i.l la'll Bv using a cylindrical lens, ue obtain a summation 
ov er i I hi; -. w e h.n c 

he i.o.i i to 


where atji is the \th element of the vector a and Mti.j) is the 
matrix element. Such a scheme for matrix-vector multiplication 
can also be used for matrix-matrix multiplication bv decompos¬ 
ing a matrix into column vectors and then multiply ing the matrix 
with each of the column vectors. 

As a result of the nature of the four-w ave mixing process in 
the medium, this matrix-vector multiplier operates on the held 
amplitudes and thus can be used to ! andie matrices and vectors 
with complex elements. It is a coherent device rather than an 
incoherent one. This aspect is distinctly different from most of 
the earlier approaches, which are all incoherent When the dev ice 
is operated in the coherent mode, the phase of each beamlct 
must be maintained uniformly over the transverse dimension of 
the beamlet. In addition, such phases must also be maintained 
fixed in the summation process. If these phases are not uniform 
over the beamlets, (he final step becomes an incoherent sum¬ 
mation as a result of the spatial averaging. Under such circum¬ 
stances. this matrix-vector multiplier operates on the intensities 
and thus handles only positive numbers 

V>e now describe a few specific examples of optical matrix- 
vector multiplication and optical matrix-matrix multiplication 
using various architectures based on tour-wave mixing and phase 
conjugation Other approaches to optical matrix processors are 
discussed in several review articles 

4.3.1. Optical matrix-vector multiplication using 
photorefractive four-wave mixing 

In the first experimental demonstration of optical matrix-vector 
multiplication using foui-wave nux-ng. phase coniugution in 
a photorefractive BaFiO. crystal w_ iseJ to perform pixel-bv 
pixel multiplication The summation required to obtain matrix 
vectoi products was performed subsequentlv by a cylindrical 
lens external to the crystal 

A schematic diagram of the experimental contiguratioi, and 
the result is shown m Fig d Mask M imprints the matrix 
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into the two counterpropagating pumping beams A collimated 
uniform beam is used as the probe. The phase-conjugate output 
( 31 , at the focal plane of the summing cylindrical lens is shown in 
the lower right. The intensity of the output spatial pattern is 
proportional to the product vector v. given by 


and mask M; imprints the vector (repeated on the mask instead 
of using a cylindrical lens to fan out the vector) 
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Fig 9 Experimental configuration and the result of optical matrix- 
vector multiplication using photorefractive four-wave mixing 
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With a slight modification of the experimental geometry. it 
has been successfully demonstrated " that both the pixel-by-pixel 
product and the summation can be done inside the nonlinear 
medium to achieve the desired matrix-vector product without 
use of an external cylindrical lens. 

The experimental geometry is illustrated schematically in 
Fig 10 The readout vector, which consists of a vertical column 
of beamlets of equal intensity, enters the crxstal at an oblique 
angle such that each beamlet traverses all of the counterpropa¬ 
gating pumping beamlets in the same horizontal plane inside the 
crystal. Since the proper elements of the matrix and vector are 
encoded in the two counterpropagating beams v ia the appropriate 
masks, the phase-conjugate output of each beamlet consists of 
the sum of the prouuct resulting from each individual encounter 
An experimental result for multiplication of a 4 > 4 binary 
matrix and a 4 x ] binary vector is shown in Fig 11 The mag¬ 
nified images of the masks used to encode the matrix, vector, 
and probe are shown in Figs. Ilia) through 1 lto. rcspectivelv 
The image representing the rinal matrix-vector product is shown 
in Fie. 1 lkii 


4.3.2. Matrix-matrix multiplication by means of color 
multiplexing 

1; ha'- been experimentally demonstrated that matrix-vector mul¬ 
tiplication using photorefractive four-wave mixing can be ex¬ 
tended to matrix-matrix multiplication h' color multiplexing 1 
The basic idea is to decompose the problem into matrix-vector 
multiplications and Curry out all of these multiplications simul¬ 
taneously in parallel by using color multiplexing This idea of 
encoding the component vectors with different colors is illus¬ 
trated in Fig. 12 for the case of 4 > 4 matrices Each row vector 
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Fig 10 Modified configuration for optical matrix-vector multiplication using photorefractive four-wave mixing In this approach, both 
the multiplication and summation are carried out inside the nonlinear crystal 
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Fig 11. Experimental results of optical matrix-vector multiplication using the approach illustrated in Fig 10. 



Fig 12 Schematic diagram illustrating the basic idea of optical ma¬ 
trix-matrix multiplication using color multiplexing 


of the matrix M> o illuminated by one color, and all of the color 
components are then combined by the prism :angular multiplex¬ 
ing into a r incie row prior to further expansion by anamorphic 
optics mot shown) to match the mask representing the second 
matrix After proper element-by-element multiplications and 
summation i summing optics omitted) into the column, the re¬ 
sulting multicolor output is demultiplexed into different color 
components that together represent the tmal product. 

4.3.3. Matrix-matrix multiplication using spatial 
convolution 

In addition to the above scheme, another method tor matrix 
matrix multiplication that uses spatial convolution by means oi 
four-wave mix me has been successfully demonstrated ~ The 
ke\ difference between this approach and those described above 
is that the nonlinear crystal is Kvatcd at the common Fourier 
plane, rather than at the common image plane, of the input matrix 
masks The encoding scheme is also ditterent from those used 
m the earlier approaches Matrix-matrix multiplication in full 
parallelism is achieved by space-multiplexing by mentis of spa¬ 
tial convolution usinc degenerate tour wave mixing 

White and Yam* have demonstrated that (spatial i convo¬ 
lution and correlation of two (iwo-dimensionali patterns can be 
achieved in real ume bv tour-wave mixing in the common Fou 
ner plane of the input patterns and recording the phase-coniugated 
output at the correspondin'. - object plane A typical experimental 
configuration •' shown scbenuticalh in fug I s I) 1 5 is a small 
• i>o;iure simulating a point source delta function), a ratten: 
represent in c -n e convolution of f t an.i I ‘2 < l I * l 2 1 i observed 
ip the output pijir. - The design o! I I anil 12 to represent -wo 
matri.e- s> . t r..,r their p-rodu. - o represented by [ |*t 2 is now 
explained 


Fig 13. Schematic diagram illustrating the basic idea of four-wave 
mixing in a nonlinear medium located at the common T ourier plane 
of the input spatial pattern. 


For simplicity. consider the simplest case of two matrices A 
and B w ith dimensions 2 x 2 and the product = AB given 
below. 



i ?41 



( 35 ) 



a l; b" * a : ;b:i anb : . * a.b; 

a ; ,bi i •* a;:b;i a ; >h ; - J;;b : ; 


(3b i 


Transparent L I c Hoist- of tour small apef'ires. each w ith ar 
intensitv transmittance proportional to each matrix element a. 

|see Fig I it o 1 1 A similar transparency 12 concspondmg tv' 
the transpose ol B tie., rows and column interchanged 1 is 
shown m Tig 14(b) Noie that the vertical separation in l 1 is 
considerablv larger than that in 12. while the horizontal elen.cnt 
separations are identical in L 1 and 12 The two-dimensional 
spatial convolution of I I and L’2 is shown in Fig 14ui. It 
consists of a total of In spots with the intensity ot each pro¬ 
portional to the cross product ot all ot the elements in A and B 
[he tour doublets Ivinc along the y -axis in Fig 14u i are drawn 
shghtlv otlset from the axis to expose the mdiv idual component 
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C = AB 

C 11 = 8 11 b 11 ' 0 12 b 21 

C 12 = a 11 b 12 * 8 12 b 22 

C 21 = a 2 1 b 11 ' e 22 b 2^ 

a 11 a " 2 ( 8 ' / b 11 b 12 . C 2 2 = *21 b 12 * 8 22 b 22 

a 2* a 22 \ b 21 b 22 ' 


Y 

* 



a (b, (ci 


Fig 14 Schematic diagram illustrating optical matrix-matrix multiplication by convolution for 2 x 2 matrices (a) Matrix A and mask U1 to 
encode A, (b) matrix B and mask U2 to encode B T , and (c) spatial pattern resulting from convolution of U1 and U2 
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Fig 15 Experimental results for optical matrix-matrix multiplication by convolution by means of four wave mixing in the spatial frequency 
domain lai and (bi are the matrix masks used for encoding the matrices A and B 1 , (c) shows the result of their spatial convolution, and (d) is 
an overexposed version of Ic' to reveal the cross terms that do not contribute to the matrix product 


In practice. the -c t v*. < > comp, •nent'- spatially oveiiap on the v- 
u\is The intensities the Mur doublets, troth top to bottom 
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Another scheme for optical matrix-vector (and matrix-matrix I 
multiplications lhai uses a phase conjugator (with a finite storage 
timei :n conjunction with a spatial light modulator iSLMi to 
eliminate the pixel-hv-pixel alignment requirement at the cost 
ot some leduction in paral'c.ism iN" 2 instead of N~) has been 
repined 0 Phase aberration due to imperfection in the optics is 
sell-corrected bv the phase-conjugation process. The optical sys¬ 
tem involved is relatively simple cor,. • ed with the other ap¬ 
proaches Figure lb is a schematic diagram illustrating the basic 
concept An SLM is used to impress the matrix and vector 
information in sequence on an input laser beam This beam is 
directed toward a self-pumped PCM, which has a finite time 
response (e g . BjTiO.j A cylindrical lens is used in the phase- 
comugutc output beam path to pertorm the summation 

The principle of operation is as follows The SLM first im¬ 
presses the matrix information onto the input laser beam This 
beam is then incident into a self pump’d PCM. which stores the 
m.ifiiv information after a unite grating formation time When 
iln matrix inlormaiior, i- removed from the SI M and all ot the 


1 ; >' M p:vel- are turned into maximum transmission condition, the phase 

1 • a civer. a'Miutj’,' beam that contain the reconstruction of the matrix 

a 'I-!-.. .1 :: ' • ml •nn.it: •; exists |.<t a t’.iute dutalior.. wltich depends on the 

. s;;; i,t: inpi' 'leas! beam, lb:;'.: th> time . if the next 

■ • ,.i live si . irtie ifn ve. lor it:! am.parallel nul 

..if' •• I'pr'. t" • ! as the piia'c s•••,!::ea‘. I'Cam plopag. 

• S; V Hems fl.v v.-.t, : :s refteser.ted as a two 
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Fig 16 Schematic diagram illustrating the concept of optical ma¬ 
trix-vector multiplication using a photorei. active phase conjugator 
in conjunction with an SLM 

dimension of the vector. A cylindrical lens in the output port 
perlorms the summation. The dark storage time during which 
the matnx information can he retrieved is determined by the 
photorefractive material and the pumping configuration. It ranges 
from seconds to microseconds, in this architecture, the same 
input beam is used for alternately writing and reading the holo¬ 
gram. 

This system can ako perform matrix-matrix multiplication by 
time multiplexing In this case, each column sector V, li = 1 
to N ’ constituting the second matrix M; is sequentially impressed 
onto the beam K- multiply with the first matrix M| according to 
the matrix-vector multiplication scheme described above To 
avoid degradation os the information of M; stored in the photore 
fractive hologram during the readout, it is necessary to refresh 
the holographic memory with M. to restore its diffraction ef¬ 
ficiency This can be d>me by reimpressing M, onto the beam 
after each readout cycle Consequently, a tola’ of 2N clock 
cycles, consisting o' N cycle"' of write and N cycles of read, 
will be required to carry out trie multiplication of two N ■ N 
matrices 

4.4. Optical interconnection 

Optical interconnection linking laser arrays and detector array" 
plays a key role in optica! computing ' Conceptually, such 
an interconnection can he achieved by using optical matrix 
\c\!e: multiplu atiou. 

i Mi i 3" i 

where v is the input vector representing the signals earned by 
an array of N lasers and v is the output vector representing the 
signals received bv the array of N detectors M is an N > N 
bmarv matrix representing the interconnection patten 

When an SLM i" used a" the matrix mask, the .nterconnectmn 
is re -.onf. curable Suppose the interconnection mask consists o! 
an \ / N array . f tran-mi"'ion v .ndows The input vector la 
row of N elementsi is tanned out using astigmatic optics so that 
the hen! Iron; each element is broadcast over a correspondin'.: 
loluinn o! window e After transmission through the intercon¬ 
nection ma-l . similar optics are used to collect the light trom 
ea.t r>w of the window > and to sum the output into a column 
vector ot \ element- Sir Is an architecture provides the nec 
esc.irv N p.iraiic!i"m However a larce traction ot energv n 
jb-orVd R . to-. ; rat.sp. em. c or Si.M W lien used as a 
swtf.r. sir. t ,ir ar. hit. . ’c-. has an er.ercv ft tic lens x "I on'. 

I k whe k r ■ t;.-.- ,i;-i . -.-ri - a-■: >! Ih.- arrav I in- iscurs b-.-c'iiise 
I k I - k of It. : • •' : er -r e . 1' e.u r- civil .cm ■ >! ti'e l rtf'at 


vector does not pass through the crossbar mask The energy loss 
increases as the dimension of the array increases and is reterred 
to as the "fan-out energy loss." For a 1000 x 1000 crossbar 
sw itch, the loss due to fan-out can be as high as 99.9 r /<. This 
is not acceptable tn high speed computing because signals are 
passing through the SLM at gigabit rates and the energy loss 
can be enormous In addition to the inherent fan-out energy loss, 
all SLMs have finite insertion loss due to imperfect transmissions 
such as those that occur from absorption and scattering, includ¬ 
ing such insertion losses, the energy efficiency to: a crossbar 
switch would become LN. where t is the transmittance (t < 1) 
of each of the windows. 

Recently, a scheme has been developed* 0 '* 1 to form a dy¬ 
namic hologram inside a photorefractive crystal such that the 
hologram will be capable of performing reconfigurable optical 
interconnection. Such a scheme utilizes nonreciprocal energy 
transfer in two-wave mixing to achieve an extremely high energy 
efficiency 

Referring to Fig. 17. we describe this new method ot recon¬ 
figurable optical interconnection. Figure 17iai describes a one- 
dimensional case to explain the concept A small fraction of the 
incident beam is split off using beamsplitter BS Tne reflected 
beam (probe beam) is then expanded by using a cylindrical lens 
before entering the SLM in this example, the input beam is to 
be connected to detectors b and d as prescribed by the SLM. 
The transmitted probe beam is then recombined with the main 
beam inside a photorefractive crystal As a result of nonrecip¬ 
rocal energy coupling, almost all of the energy in the mam beam 
is transferred to the probe beam, which carries the interconnec¬ 
tion pattern The result is an optical interconnection with a very 
high energy efficiency. 

Figure 17ih) shows the reconfigurable interconnection for 
laser arrays and detector arrays in this example (4 > 4 inter¬ 
connection!. laser 1 is to be connected with detectors b and c. 
laser 2 with detectors a and d. laser 3 with detectors c and d. 
and laser 4 with detectors a and c In terms ot matrix-vector 
multiplication, such an interconnection can be written as 
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where vq. v : . v .. and v 4 are signals carried bv the laser elements 
ot the array v 

A cylindrical lens is used to locus the two-dimensional array 
of beams into a vector Thus, as a result, detector a receives 
signals from lasers 2 and 4. detector b receives signals from 
laser I. detector c receives signal" from lasers 1. 3. and 4. and 
detector d receives signals from lasers 2 and 3. This concept 
can be extended to interconnect N lasers wnh M detectors, where 
N and M are two large numbers 

The two wave mixing described above mav be viewed as a 
real time holography in wh. h the recording and readout occur 
simuluneouslv inside the photoretra.tive crw.tl The beam¬ 
splitter and the SLM are used to generate the reference and obtect 
Reams to record a volume holoer.im that represents the intei 
connection pattern as prescribed bv the SLM I he energy cou¬ 
pling invoked in two-wave mixing ensure' lh.<: im -littraction 
efficiency during the readout o ,ilm-"i Ido Ibis requires a 
proper orientation ot the vrxsul so that the enetev ot the readout 
Ream r greatlv depleted I he high energy e!tu--iuv results tiom 
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Fig 17 Schematic drawings of optical interconnections using dynamic photorefractive holograms with N = 4. (a) 1 x N optical interconnection 
and IbINxN optical interconnection 


the fa't that most of the energy is carried by the readout beam, 
which Jot's not pass through the SLM but is diffracted into the 
interconnection pattern by the hologram 

The energy efficiency of such an interconnection pattern can 
be estimated for a crossbar switch as follows Let R be the 
reflectance of the beamsplitter It is legitimate to assume that 
the beamsplitter is practically lossless, and we may assume that 
the surface of the photorefractise crystal is antireflection-coated 
so that the 1-resnel reflection loss can be neglected Under these 
conditions, the tyyo beams that arrive at the photorefractise crys¬ 
tal have energies tl - R> and RtN. respectively Inside the 
crystal, these two beano undergo photorefractise coupling. As 
a result, most of the energy o! the pump beam H - Ri is 
transferred to the probe beam Rt \. which contains the inter¬ 
connection pattern The energy efficiency can be easily denseJ 
and is given by v 
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Fig 18 Schematic diagram illustrating the relative orientation of the 
pump beam, probe beam, and c-axis of a 30 -cut BaTi0 3 crystal for 
the measurement of energy efficiency The depleted (transmitted) 
pump beam is omitted for clarity 
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L is the interaction length, y is the coupling constant, and a is 
the bulk absorption coellictent for photoretra.tive crystals such 
as BaTiO. and SBN. the coupling constant is scry large ti e., 
■yl - li The efficiency can be written approximately as 


.votive that tor large N. the energy efliciency is limited by the 
crystal bulk absorption exp' ol.i and o maximized by using a 
beamsplitter with a very small reflectance R ti e . R (h 
In examining the energy ef tier icy. we collimate an argon 
ion laser beam at 514 > nm into 2 mm diameter by us;- ; a lens 
ot toxa! lenclh I 2 m A beamsplitter ts used to redirect 5'< 
of the energy throueh the SI M establish. the probe beam 
I he remainder "I the energy transmits through the beamsplitter 
and soc.s'i'ii'e- the pump be.i'n It- our preliminary experimem. 
the SI M w a - repi.Ked yx nti a neutral deroit. inter w ith a y artah'e 
optua; den - ;t . i- -imula'e th. tan on: energy loss [tie tw-. 
bca::_■ • k a | :< » . re-la: xx ill. ai. interaction leneti. 


of approximately 5 mm (sec Fig. ISi. The intensity ot the pump 
beam is approximately 3.2 W cm', whereas the intensity of the 
probe beam varies from 0.16 Wcm" tex 0.16 m V> cnr. The 
amplified output intensity is monitored and the result is shown 
m Fig IV. where the energy efficiency (defined as the ratio of 
the amplified output intensity to the total input intensity i is plot¬ 
ted as a function of 1 N (where N is the number ot tan-out 
channels) The solid line represents the theoretical hi using Ewp 
1 30 1 with oL and 7 L as the fitting parameters The results 
tuL = 9.2 and vL = 1 . 11 agree w ith those obtained from other 
independent measurements to within ION. Note that the energy 
efficiency is better than 24N for N = 1 through 100. In other 
words, if the lan-out loss is 99 r r (e g., a 100 > Ilk- crosshari. 
the energy efficiency of this new scheme can be 24 times better 
than that of the direct approach Bulk absorption in this particular 
crxstal accounts for approximately 650 of the energy loss 
Even though the reconfiguration time is limited by the photo- 
refractive response time, whieh is typically on the order of milli¬ 
seconds at modest intensities, the photorefractive interconnec¬ 
tion sy stem can accept very high data rate signals To demonstrate 
this tact. we frequency modulate a laser beam targon. 
\ - 514.5 nmi using an acousto-optic device to simulate a sig¬ 

nal that i' to be interconnected with some output The signal 
used had a pulse train of frequency h - 0 k>33 MHz with each 
pulse being approximately d 2 gis yy ide I his rate is clearly much 
hi ther than the reciprocal >! the photoretractive response tune 
Ilk modulated laser beam xx.o then spin into two beam'- and 
mixed in the erx stai a- described before, and the amplilicd probe 
beam was monitored witti a plioiodetedo: I he is. Bioscope 
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Fig 19. The energy efficiency as a function of 1/N, where N is the 
equivalent number of fan out channels Open circles represent the 
experimental results, and the solid line is a theoretical fit [using 
Eq (39)1 with aL and yL as the fitting parameters. 



Fig 20 (a) Signal carried by the probe beam and (b) amplified signal 
carried by the probe beam after undergoing two-wave mixing 


trace in Fig. 20u) shows the input probe signal, and the trace 
in Fig 20(h) show, the amplitied probe signal. The results show 
steads-state response in which the temporalis modulated pump 
and probe beams interact simply by diffracting off the stationary 
index grating that is created in the cry stal alter the beams interact 
for a time on the order of the photorefractise response time 
This experiment was performed morels to demonstrate the high 
signal bandwidth o! the system Optimization of the parameters 
was not done, so the result' shown in Fig 2() can cleans be 
linpri is ed 

5. COM IT SION 

We have des«.rihe.i '"me oi the interesting and important phe¬ 
nomena in pi,..; Tctra-. use rued..-. F’hotoretractive materials suOi 
a- h.LIiO SON. and BV ) are bs tar the most etlicicnt media 


for the generation of holograms and phase-conjugate waves using 
relatively low optical intensities. In addition, two-wave mixing 
in these media exhibits nonreciprocal energy transfer As a result 
of the strong energy coupling, several unique photorefractive 
phenomena occur in these media These include self-pumped 
phase conjugation, nonreciprocal energy transfer without phase 
cross talk, phase-conjugate interferometry, etc These unique 
phenomena play an important role in many of the applications 
in optical computing including matnx-vector multiplication, par¬ 
allel subtraction, reeonhgurable interconnections, etc. Some of 
these applications were presented and discussed. 
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ABSTRACT 

We propose and demonstrate the use of Fourier transform to 
achieve maximum energy efficiency in a photorefi active 
optical interconnection. The results of experimental 
investigations on reconfigurable optical interconnections 
using photorefractive holograms in a barium titanate crystal 
are presented and discussed. High energy efficiency is 
achieved by matched amplification at the Fourier plane. 
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1. Introduction 


Optical interconnection between VLSI circuits, 
computing chips or boards plays an important role in optical 
computing [1,2]. Such a scheme of computing provides the 
potential of achieving extremely high speed because it 
utilizes the fast switching of electronics and the wide 
bandwidth of optics for communication [3]. A generalized 
optical crossbar switch can provide a reconfigurable 
arbitrary interconnection (4,51 between an array of N lasers 
and an array of N detectors. Such a crossbar can be 
implemented using parallel matrix-vector multiplication 






which provides N- parallelism (5-8J. 


It is known that the N- parallelism is accompanied with 
an intrinsic fan-out energy loss [4,9], The fan-out leads to an 
energy efficiency of 1/N. We recently proposed and 
demonstrated a novel concept of reconfigurable optical 
interconnection [9-11] which can provide high energy 
efficiency. Using the beam coupling in photorefracti ve 
crystals [12.15], such a scheme of interconnection is capable 
of minimizing the fan-out energy loss and thus achieves 
extremely high energy efficiency. Although the basic concept 
has been validated and some of the preliminary results have 
been reported [9.111. the issue of energy efficiency remains a 
subject of experimental investigation. This paper reports the 
results of an extensive investigation on the energy efficiency 
of reconfigurable optical interconnection using 
photorefracti ve holograms. 

In what follows, we will briefly review the basic- 
principle of operation of the photorefractive interconnection. 
We will point out the use of Fourier transform to achieve 
maximum beam overlap in the photorefractive medium. We 
then discuss the results of our experimental investigations. 




2. Principle of Operation 


Figure 1 shows the basic idea of reconfigurable 
interconnection using photorefractive holograms. A beam 
splitter is used to split off a small portion of the beam which 
consists of beamlets each carrying information from one of 
the sources. This small portion of the beam will be called the 
signal beam. The beam splitter allows the majority of the 
energy to pass through. This major portion of the beam A 
called the pump beam. The signal beam is fanned by the 
cylindrical lens and then passes through the spatial light 
modulatorfSLM) which contains the interconnection pattern. 
After passing through the SLM. the signal beam is imprinted, 
spatially with the interconnection pattern. The pump beam 
does not pass through the spatial light modulator and thus 

suffers no energs loss due to the fan-out and all the other 

loss mechanwns associated with the SLM. 

The signal beam which contains the interconnection 
information is then recombined with the pump beam in a 
photorefractive crystal. Under the appropriate conditions, the 

signal beam can be amplified at the expense of the pump 
beam. Most of the energy of the pump beam will be 
transferred to the signal beam provided the length of 




interaction is large enough [14,15]. It is known that 
photorefracti ve two-wave mixing exhibits energy transfer 
without phase crosstalk [16]. Thus the amplified signal beam 
will also bear the interconnection pattern. The net result is 
an optical interconnection with high energy efficiency. 

Although photorefractive two-wave mixing is involved, 
such a scheme can also be considered as a real-time 
holographic interconnection with a very high diffraction 
efficiency. The combination of the beam splitter and the SLM 
in conjunction with the pump beam is used to record a 
hologram inside the photorefractive crystal. Such a hologram 
contains the interconnection as prescribed by the SLM. Once 
the hologram is formed, the pump beam can be diffracted off 
the hologram and be redirected into the arras of detectors. 
The advantage here is that the recording and the readout 
occur simultaneously. This offers the possibility of 
reconfigurab’e interconnection by 'ontrolling the SLM. 

To achieve maximum energy efficiency, it is important 
that the signal beam and the pump beam overlap completely. 
Specifically, beamlet 1 of the pump beam in Fig. 1 must 
overlap completely with column 1 of the signal beam. And 
similarly, beamlets 2. 3. and 4 must also overlap complete!) 
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with their corresponding columns in the signal beam. It is 
important to note that although the beamlets in the signal 
beam and the pump beam are intrinsically different because 
of the interconnection pattern, the individual pixels can have 
identical shape (e.g., square). Complete overlap is possible in 
the Fourier domain, provided all the pixels are identical. Let 
s(x.v) be the aperture of an individual pixel, and o( u.v ) be 
the Fourier transform. Fourier transformation of a beamlet of 
the pump beam and the corresponding column of beamlets of 
the signal beam can be performed by a lens (see Fig.2). The 
resulting amplitude distributions at the rear focai plane are 
given b; 


Pump: Pi u ,\ ) = c, u.v ) —.— (1) 

Signal: S( u.v ) = Z oi u.v ) exp 1 i o n (u.v) ) — (2) 

n 


respectively, where u.v are coordinates in the Fourier plane. 
The summation in equation (2) is over all the windows open 

in the column, and o n (u.\ ) is a phase which depends <n the 


position of the window. 








We note that the Fourier transform of a column of the 


signal beam consists of a linear superposition of identical 
patterns each with a different phase factor. Such phase 

factors exp(io n ) are due to the relative positions of the 
windows in each column. Although these phase factors ma\ 
lead to interference structures, all the energies are confined 
under the same envelop g(u,\). As a result of the shift 
invariance in Fourier transformation, each of tne signal 
beamlets overlaps completely with the pump beam at the 
F-ourier plane. Thus maximum energy coupling is achieved. 

iv_F yperi men tal I nvest; can mis 

The energy efficiency of p h o t ore f rac it v e 
interconnection depends on several parameters. In addition 
to the beam overlap mentioned earlier, it depends on crystal 
orientation, coupling constant, length of interaction, beam 
intensity ratio, insertion loss at SIM 11 MS;, and other loss 
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mechanisms such as Fresnel reflection, beam fanning [19] 
and scattering, etc. We have carried out a series of 
experiments to investigate the energy efficiency of 
photorefracti ve interconnection. These experiments are 
designed to study the dependence on those physical 
parameters mentioned above. 

Referring to Fig. 3. we consider the efficiency of 
energy transfer in photorefracti ve two-beam coupling. An 
input beam of power Pj is split by a beam sp'itter BS into a 

pump beam with power and a signal beam with power P s 
( P_ N = P, - P ). These two beams are recombined in a 
photorefractive BaTiO; crystal. The crystal is oriented such 
that the signal beam is amplified. Let P 0 be the power of the 
amplified signal beam. The energy efficiency r\ o f 

photorefractive two-beam coupling is defined as r\ - P ( , / P, . 
A theoretical expression of such an energy efficienc} was 







derived earlier in reference 9. In our experiment, we 
examine the dependence of rj on the pump-to-signal intensity 
ratio m, defined as m = P p / P s . Such a parameter can be 
changed by using a variable beam splitter. 

Using argon ion laser as the source, and a 5 x 7 x 7 
mm' BaTiO; crystal as the photorefractive medium, we have 

measured the energy efficiency (rj ) as a function of the 
intensity ratio m. The maximum efficiency is approximate!} 
10 U when the beam intensity ratio is m = 20. In addition to 
these measurements, we also repeat the experiment b\ 
inserting a neutral density filter (NDF) in the optical path of 
the signal beam. This neutral density filter is used to 
simulate the fanout energy loss. The idea here is to show that 
the amplified signal beam power does not decrease 
significantly as the input signal beam is attenuated. The 
experimental results are shown in Fig.4. We notice that an 
energy efficiency of 1.5 U is obtained when the simulated 







fanout loss is 99.99c. This is an improvement in energy 
efficiency by more than a factor of 15. As discussed earlier 
in reference 9. the optimum beam splitting ratio m (i.e. the 
pump-to-signai intensity ratio measured immediately after 
the beam splitter, and before the signal beam is attenuated 
by the NDF) depends, in general. on the photorefractive 
coupling constant yL and the fanout loss of the signal beam. 
In this particular experiment. the crystal is oriented such 
that both beams enter the crystal symmetrically at an angle 
of incidence of 20° as shown in the inset in Figure 4. This 
configuration is. however, not optimized for maximum 
energy efficiency because it only provides a coupling 
constant yL of approximately 4.5. 

The coupling constant yL depends on the crystal 
orientation relative to the beams. In another experiment, 
the same crystal is oriented such (hat the signal beam enter 
th; crystal at an angle of approximately 67' which is near 


9 







Brewster's angle. This incidence configuration minimizes the 
Fresnel reflection loss. The pump beam enters the crystal at 
an angle of incidence of 32° such that the grating induced in 
the photorefractive crystal provides a coupling constant yL of 
approximately 9.2. Again, an argon ion laser is used as the 
light source. The beam splitter has a reflectance of 95 c < 
which provides a beam splitting ratio of m = 1°. Under these 
conditions, the energy efficiency is improved significantly. 
Figure 5 shows the experimental results along with a 
theoretical fit (the solid line) using equation (?) in reference 
9. The energy efficiency is plotted as a function of (1/N) 
which is the transmittance of the neutral density filter. We 
recall that (1/N) represents the maximum energy efficiency 
achievable by a conventional \ N crossbar as a result of 

fanout loss. The orientation of the crystal relative to the 
incoming beam is shown in the inset of Figure 5. We notice 
that an improvement by a factor of 100 is achieved in energy 
efficiency when (1/N) is of the order of 1 x 10'-N Further 
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improvement in energy efficiency can be achieved by coating 
the crystal surfaces with anti-reflection dielectric films and 
by using crystals with higher coupling constants. 

In the experiments discussed above, both signal and 
pump beams are Gaussian beams with no spatial intensity 
pattern. In the interconnection applications, these beams are 
spatially modulated. As we mentioned earlier, the energs 
efficiency of the photorefractive interconnection as depicted 
in Fig. 1 depends on the spatial overlap of the beams. We 
now discuss our experimental investigations on the 
utilization of Fourier transform to achieve maximum beam 
overlap and thus to achieve maximum energy efficiencs. 

The experimental configuration is illustrated in Fig.6. 
The output from an argon laser (514.5nm) is spatia 1 -fi 1 tered 
and expanded to form a collimated beam with 2cm diameter. 
A variable beam splitter consisting of a half-wave plate and 
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a polarizing beam splitter cube is used to split the incoming 
beam into two, a pump beam and a signal beam. Another 
half-wave plate is used to rotate the polarization of the 
reflected (signal) beam by 90° from the direction 
perpendicular to the plane of the paper to the in plane 
direction. A polarizer placed downstream further filters out 
the residual perpendicular component. A 10x10 binary 
matrix mask is used to encode a spatial pattern omo the 
signal beam. The objective here is to efficiently transfer 
energy, via photorefractive two-wave mixing, from the pump 
beam to the signal beam which is encoded with the 
interconnection pattern. Any one or more (up to 100) of the 
10x10 channels can be selected by using an appropriate 
mask. 

To maximize the spatial overlap of the pump and the 
signal hp-ams inside the crystal via the properties of Fourier 
transform [see equations (lj and (2)j described in Section 2. 
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a pump mask with a single aperture identical to the unit cell 
of the signal mask is used. The two masks are placed at the 

front focal planes of two Fourier transform lenses FLj and 
FL-,. (of identical focal length f = 25cm) respectively. Here we 

use two lenses of identical focal length because optimum 
beam coupling occurs at an angle of 40“ between the beams 
which are beyound the numerical aperture of a single lenv 
The crystal is located at the common Fourier plane of the two 
masks (see Fig.6). The shift-invariant property of Fourier 
transformation ensures that, apart from a phase factor, the 
diffraction pattern from each signal channel overlaps with 
that from the pump inside the crystal. '»he intensity patterns 
of the pump and the signal at the image plane and the 
Fourier plane are shown in Fig.7. Note that the two intensity 
distributions at the Fourier plane differ significantly in fine 
structures due to multiple beam interference among the 
various signal channels. A magnified version that reveals 
these fine structures can be found in Reference [20j. 
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The key questions we address in this section are the 
following: (1) How efficiently can we transfer energy from a 
single channel to multiple channels ? (2) How does the 
energy efficiency depend on the number of signal channels 0 
and (3) How is the energy distributed among the signal 
channels ? 

To investigate how the energy efficient's depends on 
the number of signal channels, we used a variable 
rectangular aperture in front of the signal mask to vary the 
number of channels from 100 to 1. In each case, we also 
adjusted the laser power and the variable beam splitter so 
that both the pump power and the total signal input power 

are fixed at 600 p\V and 6gW, respectively, regardless of the 
number of signal channels. Figure 8 is a plot of energy 
efficiency vs. number of signal channels. The energs 
efficiency turns out to be fairly insensitise to the number (if 
signal channels. These data were taken with the geometry 








shown in the inset of Figure 4. As mentioned before, this 
particular geometry is far from optimum, and one can expect 
an energy efficiency of 30 to 409 x or higher once optimized. 

With all the 10x10 signal channels on, the intensity 
pattern at the output (imaget plane is shown in Figure 9. The 
original output pattern with the photorefractive crystal 
remo\ed is shown in Fig.9 (a). When the crystal is in place, 
the output pattern is slightly degraded as shown in Fig.9 (hi. 
The degradation is due mainh to beam-fanning [171. This can 
can be verified by observing the output patten, using a 
signal oeam which is ordinarily poianzed to-polarized i u; 
much lower in power. The results are shown in Fig.9 (c) and 
id). In botn cases the original intensity nauern is restored as 
the beam fanning effect diminished. In Fig.9 (c). significant 
reduction of beam fanning is a consequence of the much 
smaller electro-optic coefficients associate with the 
o-poLnzed beam. In Fig.9 id), it is a consequence of longer 
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time required for the fanning grating to build up when the 
beam intensity is weaker. The bottom picture in Fig.9(e) 
shows the amplified signal output when the pump beam is 
turned on. The presence of the pump beam not on!} 
increases the output energy but also significant!) reduces the 
signal degradation. To monitor the cnergv distribution 
quantitative!}, a linear detector arrav is placed along one 
rose of the signal channels at the image plane 1 he 
oscillograms representing the energx distribution among the 
10 channels are given in Fig. 10 tor the following three cases: 
(a i with the crystal removed: (hi with the crystal in phue 
and the pump beam off; tci with the crystal in place and the 
pump beam on. In all cases, the Gaussian top-hat' profile 
originated from the expanded laser beam is more or less 
preserved. and the energy distribution is sufficiently 
uniform. Note that the vertical scale of the oscillogram in 
Fig. 10(b) is 1 ()m\/division versus the 1 OOmY/dtvision scale 
in big. 10(c). The signal gain, in this case is approximate!} 20. 







Using an experimental configuration similar to the one 


shown in Fig. 6. with ail the Fourier transform lenses 
replaced b\ astigmatic optics (which image along the x-axis 
and Fourier transform along the y-a\is in the input arms, and 
\ue versa in the output arms), we also demonstrated a 6xb 
generalized <.rossbar-sw itch. An appropriate pump mask and 
a signal mask for a specific interconnection pattern are 
show t. in the lelt hand side ot Fig.II. The corresponding 
c tetisit) patterns inside the crwal at the Fourier plane are 
shown on the right. ldeallx. we want each pump channel to 
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4. Summars 

In summary, we have proposed and demonstrated the 
use o* Fourier transform to achieve maximum energy 
efficiency in a photorefractive optical interconnection. We 
iia' e shown that the energy efficiency of photorefractive 
dynamic holograms can be of the order of 30 f f to 4(l r ; . hi 
addition, we have also verified that the energy efficiency is 

insensitive to the number of signal channels, and that the 

energy distribution among the signal channels is sufficiently 
> ,n itorm. Finally we have described a 6x6 crossbar-switch 
using this approach:. 

All the experiments described in this paper were 
carried out with photographic masks replacing the SLMs. To 
demonstrate the reconfigurability and to studs the 

reconliguration time, one can no longer avoid the use of SI. Ms 

which, in general, have relatively poor contrast ratio (of the 






order of 10 to 30) relative to the photographic masks. We 


are current!) investigating the effect of limited contrast ratio 
of SLMs on the operation and the reconfiguration time. The 
results will be the subject of future publications. 
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FIGURE CAPTIONS 


Fig.l A schematic diagram illustrating the basic idea of a 

reconfigurable optical interconnection using a 
photorefractive crystal in conjunction with a spatial light 

modulator. 

Fig.2 A schematic diagram illustrating the basic idea of using 

Fourier transform to achieve the spatial overlap of the pump 
beam and the signal beam. 

Fig.3 An experimental configuration for measuring the energy 
efficiency in photorefractive two-beam coupling. 

Fig.4 The energy efficiency (rj ) of two-beam coupling in a 

photorefractive barium titanate crystal a^ a function of 
pump-to-signal power ratio (m). The symbols "0". "O". 

and "A”, represent the cases without any neutral densitx 
filter, with neutral density filters NDF1, NDF2. and NDF3. 
respectively, in the signal input arm. The percentage 
transmittance of the three neutral densitx filters are 12.2. 
0.74 f /i, and 0.129f, respectively. The orientation of the 
crystal relative to the beam is given in the inset. 

Fig.5 Energy efficiency (r)) of photorefractive two-beam coupling 
in a barium titanate sample as a function of the 
transmittance of a neutral density filter placed in the signal 

input arm. The transmittance is labeled 1/N to relate it to 
the fanout loss of a NxN permutation crossbar network. The 






orientation of the crystal relative to the beam is given in the 
inset. 

Fig.6 An experimental configuration for a 1-to-NxN (for N = 10) 
broadcasting network using photorefractive holograms at 
the Fourier domain. 

Fig.7 The intensity patterns of the masks for the probe and the 

pump beams at the image plane and the Fourier plane. 

Fig.8 Energy efficiency (r] ) as a function of number of signal 

channels (N) in a 1-to-NxN broadcasting configuration using 
photorefractive holograms in a barium titanate sample. 

Fig.9 A 10x10 intensity pattern at the output image plane under 

various conditions: (a) with the crystal removed, (b) with 
the crystal in place and pump beam off. total power = 

400|iW, e-polarization, (c) same as in (b) but with 
o-polarization, (d) same as in (b) but with total power 
= 3p\V, (e) amplified output signal, total signal input = 3u\V. 
pump input = 600gW . 

Fig. 10 The intensity distribution of a selected row of 10 out of the 
10x10 channels shown in Fig.9. (at with the crystal 
removed, (b) with the crystal in place and pump beam off . 
total input power =10p\V, (c) with the crystal in place and 
pump beam on , input pump pow'er = 400pW. The vertical 
scale (per division) in the pictures are 50mV. lOmV. and 
lOOmY, respectively. 





Fig. 11 The intensity patterns of the masks for the probe and the 
pump beams at the image plane and the crystal plane for a 
6x6 generalized crossbar network. 
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On Image Amplification by Two Wave Mixing 
in Photorefractive Crystals 


John H. Hong, Arthur E. Chiou, and Pochi Yeh 
Rockwell Science Center 
1049 Camino dos Rios, A25A 
Thousand Oaks, CA 91360 


Abstract 

Coherent amplification of images by two wave mixing in photorefractive crystals is 
examined with attention given to processing in the Fourier domain. It is shown that the 
gain that is experienced as the probe image traverses the crystal is uniform across the 
image. The gain can be expressed as a well known function of the average probe to pump 
intensity ratio. Experimental verification is given to support the theory. 





The phenomenon of energy’ transfer in two wave mixing in nonlinear optical media 
has been used effectively for the coherent amplification of images [1.2]. In diffusion 
driven photorefractive crystals, for example, it is known that the inherent 90° phase shift 
between the induced refractive index grating and the holographic intensity grating is 
responsible for the energy coupling. Shown in Fig. 1 is a typical image amplification 
experiment in which, with the correct geometry', the signal image wave becomes amplified 
at 'he <*xnense of the pump wave intensity. The model that is currently being used to 
describe coherent image amplification is a simple one developed for the case of mixing two 
plane waves. Since the signal image wave can be expressed as a superposition of many 
plane waves whose propagation directions vary' over a limited range (this range is small 
when compared to the difference between the average signal wave direction and the pump 
wave direction), the simple plane wave interaction model must be modified to be more 
accurately applicable. In this letter, we present a simple modification of the theory of two 
wave mixing involving image bearing beams and discuss in particular the amplification 
process when the crystal is placed at the Fourier transform plane of lens LI. Experimental 
evidence is also given to support the theory. 

Shown in Fig. 2 is the image amplification geometry in which the amplitudes of the 

image bearing waves are represented by {Ej}, j=1.2.N, and the pump plane wave is 

denoted by Ep. In photorefractive crystals, the index modulation depends on the intensity 
interference pattern so that the present analysis begins with the intensity expression which 
is given by 

1 r N n N T 

I(r) - Ip < 1 + t— Re Ep^ E j^expfLkp-krp)-rJ + ^ X EqE m exp[i(kq-k m )-r) .,(1) 

^ m=l q^m m=l 

where k p and k m are the wave vectors of the pump and the mth probe beams expressable, 
according to Fig. 2, by 






kp = ^ (-£sin0 + fcosQ ) , k m = t^( £sin8 m + £cos0 m ) , (2) 

* /v A 

N 

and the total intensity is Io = !E n l- + X IE m l~. If the probe waves span a sufficiently small 

m=l 

angular spectrum (much smaller than the angle between the pump and the average probe 
direction such that N’58 « 0), then the interference between the probe waves will induce 
index gratings which are much weaker than those written by the interference between the 
probe waves and the pump wave [3]. This is the scenario of interest, and the index 
distribution can then be approximated as 

f N 

n = n 0 - Re J i -m- X E ;1 E m exp[-j(k p -k m ) rj l (3) 

m-i 


where n \ is the index modulation coefficient which is assumed to be the same for all of the 
gratings. The rc/2 phase shift is due to the diffusion driven process. With this and the 
assumption of slowly varying amplitudes for each plane wave involved in the interaction. 
Maxwell’s equations yield the following coupled system: 


dEp 

dz 


12 . V ip i 2 p 
Iq ^ IL m' L r- 


m = 1 




ni>. 

4cos8 


dEp, 

dz 


7 ^- !EJ- E„,. m=1.2.N, 

<0 '■ 


n]X 

4cos0„, 


(4) 


The standard separation of intensity and phase equations is useful, and in particular, the 
intensities of the pump and probes, I p =lEpl- and I m =lE m l 2 . are governed by 



E P = 2y p 


r = Oy 

1 m - ;m* 


(?) 








In the assumed diffusion driven scenario, the phases of the beams are decoupled so 
that the intensity equations describe the two wave mixing process completely. To obtain 
the final desired result, we make one further approximation which is consistent with the 
assumptions made thus far. We have assumed that 0m= 6, so that it is justified to 
approximate IT, = r m . With this, the equations become 


*k = _Ly i i 

dz Io £ '=>■ 

.N- 


These equations are easily integrated to yield the following solutions. 


S In exnt-Fz) 

I (z , = --- , 

1 + P exp(-Tz) 

Fq’z) _ 1 + p 

" 1 + (J exp(-Tz) 


P=V 


X I m <0) 


The above is the main result of this letter. The solutions indicate that the probe waves are 
amplified uniformly with a gain factor which is a function of the ratio of the pump beam to 
the total intensity of the probes. Small deviations from this are expected as the coupling 
coefficient is a function of the actual probe beam angle, and also as harmonic distortions of 
the gratings result from the modulation depth approaching unity [4-6]. 

The uniformity of the amplification has been verified by the following experiment 
(see Fig. 3). An output beam from an argon laser (514.5 nm) is spatially filtered and 
expanded to pass through a rectangular aperture (RA: 2mm x 6mm). A variable beam 
splitter consisting of two half-wave plates and a polarizing beam splitter cube is used to 
split the incoming beam into two (the pump and 'he probe) and to vary the intensity ratio of 
the pump and the probe beams; the polarizer positioned after the half wave plate in the 









reflected path ensures that both transmitted and reflected beams have the same 
polarizations. The probe beam is transmitted through five rectangular windows (0.5mm x 
2mm each) in a mask (MAj) to form 5 probe beamlets. In the other arm, the pump beam 
illuminates an identical rectangular window in a second mask (MAi). Two spherical lenses 
(FL] and FL2), one in each arm. are used to Fourier transform the two spatial patterns onto 
the crystal plane. A spherical lens (FL 3 ) is used to re-image the spatial pattern carried by 
the amplified probe beam onto the detector plane where the optical intensity profile is 
sampled by a linear detector array and monitored by a storage oscilloscope. We have 
v aried the pump-to-total probe power ratio in the range 1000 ~ 1 and observed that the 
energy transferred from the pump is equally distributed among the probe beamlets to within 
20 r ;-. A typical example of the intensity profiles is shown in Fig. 4. The lower trace is the 
intensity profile of the five probe beamlets transmitted through the crystal when the pump is 
off. The upper trace represents the corresponding profile when the pump beam is turned 
on. For this specific case, the optical powers of the pump and the total probe are 0.6 m\V 
and 0.2 mW. respectively. When all but one of the probe beamlets are blocked, almost all 
of the depleted pump energy, originally distributed among the five channels, redirects itself 
into the single unblocked channel. 

The dynamics of grating adjustments made by multiple probe beams were observed 
in the following experiment. A BaTiOy crystal was placed at the intersection of two weak 
probe waves and a strong pump beam as shown in Fig. 5. As shown, the two probes were 
arranged so that they propagate with a small angle (~ 1 degree) between them in 
comparison with the angle between them and the pump beam (~ 30 degrees). The results 
are shown in Figs. 6 a-b. The input intensities of the two probes were roughly the same 
fP~10|iW) with the input pump power set at P=5mW. Shown in Fig. 6 a are the probe 
intensities (after passage through the crystal as seen by detectors 1 and 2 ) as they are 
strobed on with shutters. The amount of pump power that is scattered into the probe 
detectors was small compared to the unamplified probe intensities and correspondingly 


D 







becomes negligible when compared to the amplified probe beams. Fig. 6b shows the 
amplified probe beams as the probes were strobed by the shutters (top trace=probel, 
bottom traee=probe2; scale factor=l V/mW). The power sharing is clearly evident in the 
mid portion of the trace where both probes are on. \ v hen the order of events were changed 
so that probe 2 was turned on first, the results were the same so that no hysteresis effects 
were seen. Note that probe 1 is slightly stronger than probe 2. This is due to the fact that 
the two beam coupling constants (one for each probe) were not the same. 

In conclusion, we have analyzed the uvo wave mixing phenomenon when a 
multitude of probe beams are used. The theoretical results shown indicate that the pump 
power is shared by the probe beams. Experimental evidence was given to validate this 
power sharing effect. 

We acknowledge useful discussions with our colleagues, Fred Vachss and Ragini 
Saxena. Also, at the time of the submission of this paper, work [7) presented at the OS A 
annual meeting and another recently published |S| which dealt with similar subjects were 
brought to our attention. This work is supported primarily by the Defense Advanced 
Research Asencv. 
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Figure Captions 

Typical image amplification experiment using two wave mixing in photorefractive 
crystal. 

Two wave mixing with the crystal at the Fourier plane of the input image. The pump 
Ep is a plane wave and the image effectively consists of a set of point sources. 
Experimental setup used to verify uniform gain. MAI consists of five identical 
apertures and comprises the probe image to be amplified. Because of the gaussian 
illumination, the intensity pattern at MAI is not uniform. MA2 is the pump image 
aperture. 

Input dower amplitude) and amplified (higher amplitude) intensity distributions. Each 
portion of the overall input distribution receives roughly the same gain. 

Experimental setup used to record the dynamics of multiple grating buildup. 

Dynamics of two gratings in the same volume, a) input probe intensities (upper 
corresponds to «1 and lower -2 in reference to Fig. 5). b; output amplified intensities 
(history: -1 turned on. *1 turned off. #2 turned on. #2 turned off. #1 turned on. =2 
turned on, *1 turned off, 
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Fig.. 6 (b ) 
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Reconfigurable optical interconnection using photorefractive holograms 

Arthur Chiou, Pochi Yeh, Scott Campbell, and John Hong 

Rockwell International Science Center 
1049 Camino Dos Rios, Thousand Oaks, CA 91360 

ABSTRACT 


The use of photorefractive holograms in conjunction with a spatial light modulator (SLM) to 
realize a reconfigurable optical interconnection with very high energy efficiency is briefly 
discussed. In this approach, the SLM is used as a programmable binary matrix mask to encode 
the interconnection pattern to a coherent laser beam whereas the photorefractive crystal is 
used as a dynamic holographic medium to store the pattern and to efficiently diffract the 
readout beam into the selected channels at high speed. Vie report recent experimental results 
on issues such as energy efficiency, reconfiguration time, contrast, and uniformity. 

1. INTRODUCTION 


Optical interconnection between VLSI circuits, computing chips or boards plays an 
important role in parallel optical computing. 1 , 2 Such a scheme of computing provides the 
potential of achieving extremely high speed because it uses the fast switching of electronics and 
the wide bandwidth and massive parallelism of optics for communication. 3 A generalized 
optical crossbar switch can provide an arbitrary and reconfigurable interconnection 1 ', 5 ' between 
an array of N lasers and an array of N detectors. Conceptually, such an interconnection can be 
achieved by using optical matrix-vector multiplication, 51 - 8 

V' i MV, (1) 

where V is the input vector representing the signals carried by an array of N lasers and V' is the 
output vector representing the signals received by the array of detectors. M is an N * N binary 
matrix representing the interconnection pattern. Such an architecture provides the N 2 
parallelism required for arbitrary interconnection. When an SLM is used as the matrix mask, 
the interconnection is reconfigurable. 

The N 2 parallelism is accompanied with an intrinsic fan-out energy loss; 8 , 9 the fan-out 
ieads to an energy efficiency of 1/N. For interconnection with large N (e.g., N - 1000), this 
large amount of energy loss is intolerable. We recently proposed and demonstrated a novel 
concept of reconfigurable optical interconnection 9 ,’* 0 which can provide very high energy 
efficiency. Using the nonreciprocal energy coupling in photorefractive crystals, 11 - 19 such a 
scheme of interconnection is capable of minimizing the fan-out energy loss; thus, it achieves 
extremely high energy efficiency. Although the basic concept has been validated 9 , 10 using 
fixed binary matrix masks (or transparencies) and the energy efficiency has been measured and 
reported, 15 , 16 the demonstration of reconfigurability with high energy efficiency and high 
signal-to-noise-ratio (S/N) using an SLM remains a subject of experimental investigation. This 
paper reports the results of the investigation on the contrast of reconfigurable optical 
interconnection using a liquid crystal television (LCTV) in conjunction with a photorefractive 
barium titanate crystal. The effect of finite contrast of the SLM on the performance of the 
holographic interconnection is also discussed. 

In what follows, we will briefly review the basic principle of operation of the 
photorefractive interconnection. We will summarize some of our earlier results on energy 
efficiency measurement and on the use of Fourier transforms to achieve maximum beam 
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overlap in the photorefractive medium. 16 We then discuss the results of our experimental 
investigations on the issue of contrast. 

2. PRINCIPLE OF OPERATION 


Figure 1 shows the basic idea of reconfigurabie optical interconnection using photore- 
x 'CCti'»'£ F.jlogi uiY. j. A twain wpLUci is used to Split off a small portion of the Deam wnicn 
consists of beamlets each carrying information from one of the sources. This small portion of 
the beam will be called the signal beam. The beam splitter allows the majority of the energy to 
pass through. This major portion of the beam is called the pump beam. The signal beam is 
fanned by the cylindrical lens and then passes through the spatial light modulator (SLM) which 
contains the interconnection pattern. After passing through the SLM, the signal beam is 
imprinted spatially with the interconnection pattern. The pump beam does not pass through the 
spatial light modulator and thus suffers no energy loss due to the fan-out and all ^>ther loss 
mechanisms associated with the SLM. 



V Vv 


Fig. 1 A schematic d.ag-arn illustrating the basic idea of a reconfigurabie optical 
interconnection using a photorefractive crystal mconjunction with a spatial light modulator. 

The signal beam which contains the interconnection information is then recombined with 
the pump bean, in a photorefractive crystal. Under the appropriate conditions, the signal beam 
can be amplified at the expense of the pump beam. Most of the energy of the pump beam will be 
transferred to the signal beam provided the length of interaction is large enough. 1 V 1 * 

Although photorefractive tw'o-wave mixing is involved, such a scheme can also be 
considered as a real-time holographic interconnection with a very high diffraction efficiency. 
The combination of the beam splitter and the SLM in conjunction with the pump beam is used to 
record a hologram inside the photorefractive crystal. Such a hologram contains the 
interconnection as prescribed by the SLM. Once the hologram is formed, the pump beam can be 
diffracted off the hologram and be redirected into the array of detectors. The advantage here 
is that the recording and the readout occur simultaneously. This offers the possibility of 
reconfigurabie interconnection by controlling the SLM. 

To achieve maximum energy efficiency, it is important that the signal beam and the pump 
beam overlap completely. Specifically, beamlet 1 of the pump beam in Fig. 1 must overlap 
completely with column 1 of the signal beam. And similarly, beamlets 2, 3, and 4 must also 
? overlap completely with their corresponding columns in the signal beam. Although the beamlets 

in the signal beam and the pump beam are intrinsically different because of the interconnection 
pattern, the individual pixels can have identical shape (e.g., square). Complete overlap is 
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possible in the Fourier domain, provided all the pixels are identical. The results of our 
experimental investigations of beam overlap and the energy efficiency were presented 
recently 15 and discussed in detail in a separate paper. 16 The key results pertinent to the 
energy efficiency are summarized in the next section. 

3.0 ENERGY EFFICIENCY 


f 


# 


t 


We first investigated the energy efficiency due to beam coupling in a photorefractive 
crystal. Referring to Fig. 2, we define the efficiency (h) of energy transfer in photorefractive 
holograms as n ^ P 0 /P,, where Pj and P Q are the optical power of the input beam and the 
amplified signal beam, respectively. A theoretical expression of such an energy efficiency was 
derived earlier in reference 9 Using an uncoated barium titanate crystal (~ 5*7.7 mm, 
30°-cut, from Sanders Associates) as the photorefractive medium and an argon ion laser 
(514.5 nm) as the source, we have achieved an energy efficiency as high as 30%. Further 
improvement in energy efficiency can be achieved by coating the crystal surfaces with 
antireflection dielectric films and by using crystals with highei o,i l pi... fc euiislarus. Ultimately, 
the energy efficiency is limited by the bulk absorption in the crystal. 



INPUT BEAM 


r 


. PUMP BEAM 



AMPLIFIED 
SIGNAL OUTPUT 

* 

ENERGY EFFICIENCY (r;) . -_2 

P. 

PROBE-TO-PUMP POWER RATIO (r) « 


Fig. 2 An experimental configuration for measuring the energy efficiency in photorefractive 
two-beam coupling. 

Using an experimental configuration as schematically illustrated in Fig. 3, we have also 
successfully demonstrated 15 , 16 that : 

• the energy efficiency is insensitive to the lateral position of the signal window in the 
SLM (shift-invariant). 

• the energy efficiency is insensitive to the number of signal channels 

• the energy distribution among all the signal channels is fairly uniform. 

The properties listed above are due to matched amplification at the Fourier plane. 16 As 
indicated in Fig. 3, each individual pixel of the signal beam is identical to that of the pump 
beam to ensure maximum beam overlap at the Fourier plane. Another desirable characteristic 
of two-wave mixing at the Fourier plane is the preservation of contrast of the input pattern. In 
the next section, we discuss the issue of contrast and a simple method to enhance the contrast. 
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Fig. 3 An experimental configuration for a 1-to-N » (S' (for N = 10) broadcasting netw-ork i.K'pg 
ohotorefractive holograms at the Fourier domain. 

4.0 CONTRAST 


All of our experimental results reported earlier and summarized in the previous section are 
t obtained with photographic masks replacing the SLMs. To demonstrate the reconfigurability 

and to study the reconfiguration time, an SLM in which the interconnection pa item can be 
varied either electrically or optically can be used. In our experiment, we used * liquid crystai 
television (LCT\) which has a relatively poor contrast ratio (of the order of - to 3C)‘V 8 
compared with photographic masks. In the parallel matrix-vector multiplication approach, the 
maximum number of fan-out of the optical cross-bar can be limited bv the finite contrast of the 
SLM. 

As a simple example, consider an M * M permutation matrix mask with contrast ratio 
C - T j/To, where T ; and T c , are the intensity transmittance of the "ON" state and the "OFF" 
state, respectively. Referring to Fig. 4 (illustrated for the case of M = 4), let us assume that 
the signal from each source is intensity modulated between zero and peak value P. The signal 
level received by each detector through the "ON" cell is PT|, whereas the noise level received 
through the (M-l) "OFF" cell is P(M-l) T 0 . The signal to noise ratio S/N is therefore given by 

S/N’ = Tj/(M-I) T 0 = C/(M-i) = C/M for M » 1. (2) 

The condition S/N > 1 is satisfied provided t.iat the number of fan-out channels (M) is iess than 
the contrast ratio (C). 

* A simple method to improve the contrast is by double-passage through the same SLM via 

retroreflection or phase conjugation. In principle, if the contrast ratio for single passage is C, 
the contrast ratio for double passage will become C 2 . The experimental configuration used to 
9 verify this fact is shown schematically in Fig. 5. An expanded and collimated laser beam from 

an argon laser (514.5 nm) is transmitted through a circular aperture (diameter = 14 mm) located 
at the input plane. Half of the aperture is masked by a neutral density filter (NDF) to form a 
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2 D SPATIAL LIGHT MODULATOR 



Fig.4 Reconfigurabk interconnection using an SLM with finite contrast. 



Fig. 5 An experimental configuration for contrast-enhancement by double-passage. 

binary amplitude mask at the input plane. The transmitted beam is collected by a lens (L) and 
directed into a barium titanate crystal. The crystal is oriented so that the incoming beam is 
retro-reflected via self-pumped phase conjugation. 19 After traversing through the aperture, 
the phase conjuga'ed beam is sampled by a beam splitter through imaging optics which reimages 
the binary mask on the output plane. The input intensity distribution (i.e. single-massage 
through the mask) and the output intensity distribution (i.e., double-passage through the mask) 
are measured by scanning a detector with a small aperture (diameter *= 1mm' across the beam 
diameter at the input and the output planes, respectively. The experimental esults(for a mask 
with optical density OD = 0.5) are shown in Fig. 6(a) and (b). Apart from some imperfection at 
the edge of the output image, the output contra" ratio is approximately the square of that of 
the input as expected. Instead of a phase-conjugate mirror, a plane mirror can also be placed at 
the input plane immediately after the mask to retro-reflect the beam back through the binary 
i mask. The resulting intensity distribution is shown in Fig. 6(c). 'Je have repeated the same 

measurement with a different NDF (OD = 1.0) and also with a LCTY at the input plane; the 
results are tabulated in Table 1. Photographs of the output image of an arbitrary binary pattern 
written on a LCTY are shown in Fig. 7(a) for mirror reflection and (b) for seif-pumped phase 
conjugation. The phase conjugation has the key advantage of being intrinsically self-aligned. 
The grating formation time for seif-pumped phase conjugation using a typical barium titanate 
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Fig .b Intensity pro:nos arid contrast of a bmar> image: (a) at the input plane, (b> at the output 
plane after noun.e-passage via phase-conjugate reflection, (c) at the output plane after double¬ 
passage via mirror reflection. The arrows indicate the intensities (of the bright and dark 
regions) used for calculating the contrast ratio. 

Table ! 

Cortrast Enhancement of Binary Images by Double-Passage 

CONTRAST RATIO OF VARIOUS IMAGES 
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crystal a: a feu tens of milliwatt of optical power is, however, relatively slow ( a fraction of a 
sec o' o to several seconds). Note that the output obtained by mirror reflection is fairly uniform 
(see Fig.b(c) and Fig. 7(a)) and the reflection is practically instantaneous. 

The photorelracuve interconnection holograms used to diffract the pump Dean: into the 
signal channels specified b> the mat: ,x mask do not degrade the contrast provided that the 
pmco-efra' live crystal is planed at the Fourier plane. 1 In contrast, photorefractive image 
amplification via two-wave mixing in the image plane tends to decrease the contrast of the 
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Fig.7 Double-passage image of a binary matrix pattern or. a liquid crystal T.V. (a) via mirror 
reflection, (b) via phase-conjugate reflection. 

input image. 2 ', 1 : An expenmenta! result varifying this contrast preserving property is given in 
Fig. 8. Intensity profiles of a binary image and its amplified output are shown in (a) and (b), 
respectively. Note that the contrast ratio of the input and the output image are almost 
identical and that the vertical scale in (b) is 20 times that of (a). 

SC50017 



(a) 


(b) 


Fig.8 Experimental results showing contrast-preserving image amplification by photorefractive 
two-wave mixing at the Fourier domain; (a) input intensity profile, (b) amplified output intensity 
profile. 

To extrapolate from the above discussion, if an SLM with a contrast ratio of the order of 
100 (for example, a ferroelectric liquid crystal SLM) is used in the double passage mode, the 
maximum number of channels (M) of a M > M crossbar switch limited by the contrast ratio of 
the SLM can be as high as 1 0 1 *! 

5.0 EXPERIMENTAL DEMONSTRATION OF RECONFIGURABILITY 


Using an experimental configuration similar to the one illustrated in Fig. 3, with the 
photograpmc mask (transparency) r eplaced by an SLM, we have demonstrated the reconfigur¬ 
ability of the interconnection by writing new interconnection patterns on the LCTV. The 
patterns are generated by an IBM PC and sent to the LCTV. With optical power of the order of 
0 tens of milliwatts, the frame rate is currently limited by the photorefractive response time to a 
few frames per second. The data rate, however, is not limited by the photorefractive response 
time and can be higher than several megahertz. 9 To enhance the contrast of the binary pattern, 
9 the signal beam is passed through the SLM twice via a retroreflecting mirror (as described in 
Section 4.G) prior to interacting with the pump beam in a photorefractive barium titanate 
crystal. 
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6. SUMMARY 

We have described a novel method for reconfigurable interconnections using 
photorefractive holograms in conjunction with an SLM. Several important features, including 
high energy efficiency (> 30%), uniform energy distribution, and high data rate transmission 
(> several megahertz) of this dynamic holographic interconnection are discussed. We have also 
described the problem associated with the finite contrast of the SLM and have shown that the 
contrast can be greatly enhanced (i.e., squared) by double passage through the SLM either by a 
retro-reflecting mirror or by a phase-conjugator. We also report the demonstration of 
reconfigurability by using a LCTV in conjunction w'ith a photorefractive barium titanate crystal. 
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